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ABSTRACT 


This  research  represents  the  first  phase  of  a new  attack  on 
the  difficult  problem  of  predicting  the  results  of  a collision 
between  polyatomic  molecules.  We  treat  the  translational  motion 
of  the  molecules  by  classical  mechanics,  and  the  internal  motions 
by  quantum  mechanics.  The  effect  on  the  trajectory  of  a quantum 
transition  of  the  internal  coordinates  is  obtained  by  requiring  that 
the  classical  conservation  laws  of  energy  and  angular  momentixm 
be  satisfied  at  all  times.  This  method  gives  a description  of  the 
molecular  trajectories  in  both  adiabatic  and  non  adiabatic 
collisions  which  is  much  more  detailed  than  has  been  considered 
previously.  Such  detailed  information  about  the  trajectories  is 
especially  important  for  applications  to  the  kinetic  theory  of  gases. 

Our  theory  requires  that  one  know  the  relative  probability 
that  a quantum  transition  of  the  internal  coordinates  takes  place 
at  various  points  along  the  collision  trajectory.  Because  of 
diffraction  effects  and  limitations  due  to  the  uncertainty  principle, 
such  information  is  not  directly  available  from  quantum  mechanics. 

* This  work  was  carried  out  at  the  University  of  Wisconsin 
Naval  Research  Laboratory  and  was  supported  in  part  by 
Contract  N7onr-E851i. 


An  approximate  formula,  suggested  by  the  form  of  the  quantum 
mechanical  equations,  is  proposed  for  calculating  this  probability 
distribution. 

The  theory  is  illustrated  by  sample  calculations  for  the  exchange 
of  rotational  energy  between  linear  polar  molecules.  On  the  basis 
of  some  of  these  calculations  and  some  earlier  work  by  F.  J.  Krieger 
j~Proj.  RAND  Report,  RM-646  (195})J  a method  is  proposed  for 
calculating  the  viscosity  of  polar  gases.  The  gas  is  viewed  as  a 
mixture  of  molecular  species,  with  each  set  of  internal  quantum 
numbers  defining  a separate  species.  In  non  resonant  collisions  be- 
tween linear  molecules,  the  interaction  energy  is  taken  to  be  of  the 
form  of  the  Lennard -Jones  potential,  with  a repulsive  energy  term 
proportional  to  the  minus  twelfth  power  of  the  separation  and  an  at- 
tractive term  proportional  to  the  minus  sixth  power  of  the  separation. 
When  the  molecules  are  in  resonance,  an  additional  term  proportional 
to  the  minus  third  power  of  the  separation  must  be  included.  This 
term  may  be  either  repulsive  or  attractive;  its  coefficient  can  be 
j calculated  from  a knowledge  of  the  internal  quantum  numbers. 

I 

I In  order  to  test  this  proposed  method  for  predicting  the  viscosity 

I of  polar  gases  it  would  be  necessary  to  calculate  the  dynamics  of 

! collisions  governed  by  the  resonant  potential  in  which  ihe  term  in 

the  cube  of  the  separation  corresponds  to  a repulsive  force.  Calcu- 
lations for  the  Lennard-Jones  potential  and  for  the  resonant  potential 
in  which  the  cubic  term  corresponds  to  an  attractive  force  are 
■ already  in  existence. 
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I.  BACKGROUND  FOR  THE  STUDY  OF  INELASTIC 
MOLECULAR  COLLISIONS 


1 - 1 Introduction 

The  kinetic  theory  of  dilute  gases,  as  developed  by  Boltzirann, 
Enakog,  Chapman\  and  others,  is  valid  only  for  the  case  of  mon- 
atomic gases,  since  the  potential  energy  of  interaction  between  two 
molecules  is  assumed  to  depend  only  on  the  distance  of  separation. 
The  theory  of  transport  phenomena  in  gases  made  up  of  molecules 

with  internal  degrees  of  freedom  has  recently  been  developed  by 

2 3 

Wang  Chang  and  Uhlenbeck  , and  by  J.  de  Boer  . 

Whereas  the  transport  coefficients  in  the  Enskog  theory  may  be 

calculated  from  a knowledge  of  the  angle  of  deflection  in  a bimolecu- 

lar  collision,  the  theory  for  polyatomic  molecules  requires  knowledge 

of  the  full  set  of  differential  scattering  cross  sections  corresponding 

to  all  possible  quantum  transitions  of  the  internal  degrees  of  freedom. 


1.  The  classic  text  on  the  subject  is:  S.  Chapman  and  T.  G. 
Cowling,  The  Mathematical  Theory  of  Non  Uniform  Gases, 
Cambridge  University  Press,  1939. 

A more  recent  text  o£  wider  scope  which  gives  up  to  date 
devebpments  in  the  field  and  emphasises  the  practical 
application  of  the  theory  is;  J.  O.  Hirschfelder,  C.  F.  Curtiss, 
andR.  B.  Bird,  The  Molecular  Theory  of  Gases  and  Liquids, 
Wiley,  1954.  This  reference  will  be  quoted  frequently  with 
page  references,  and  will  be  abbreviated  MTGL. 

2.  C.  S.  Wang  Chang  and  G.  C.  Uhlenbeck,  Univ.  of  Michigan, 
publication  CM-681  (1951). 

3.  Unpublished,  See  MTGL.  p.  501. 
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Unfortunately,  it  has  not  been  possible  to  evaluate  these  differential 
scattering  cross  sections  for  a single  case,  because  of  the  complexity 
of  the  quantum  naechanical  description  of  a thermal  collision  between 
molecular  systems.  The  work  reported  in  this  thesis  is  directed  to- 
ward an  evaluation  of  these  differential  scattering  cross  sections  by 
means  of  a semi-classical  description  of  the  collision  process. 

The  semi-claaaical  formulation  o‘.‘  the  collision  problem  de- 
pends on  the  assumption  that  the  relative  translational  motion  of  the 
colliding  molecules  may  be  described  by  classical  mechanics,  while 
the  internal  motions  - electronic,  vibrational,  and  rotational  - are 
described  by  quantum  mechanics.  This  method,  sometimes  called 
the  method  of  impact  parameters,  has  frequently  been  applied  to 
cjllision  problems  in  which  it  was  desired  to  determine  only  the  total 
probability  of  an  internal  transition,  and  in  which  the  classical  tra- 
jectories couid  be  approximated  by  straight  lines  along  which  the 
molecules  moved  at  constant  velocity. 

In  this  research  we  describe  an  extension  of  the  semi-classical 
method  which  gives  the  detailed  information  about  the  trajectories, 
which,  together  with  the  probability  of  internad  transitions,  is  neces- 
sary for  evaluation  of  the  differential  scattering  cross  sections.  The 
extension  depends  on  the  assumption  that  the  collision  is  nearly 
adiabatic,  so  that  the  Schroedinger  equation  describing  the  internal 
motions  may  be  solved  for  every  inter  molecular  separation  with  that 
quantity  appearing  only  as  a parameter  in  the  wave  functions  and 
energies  of  the  solutions. 

The  energy  term  so  calculated,  depending  on  the  intermolecular 
separation  and  internal  quantum  states,  is  taken  as  the  effective 
potential  energy  of  interaction  from  which  the  classical  trajectories 
may  be  calculated  by  integration  of  the  equations  of  motion.  In 


principle  chis  integration  can  always  be  carried  out,  though  in  most 
cases  it  will  be  necessary  to  resort  to  numerical  methods.  The 
"elastic"  trajectories  so  calculated  describe  the  motion  of  the  mole- 
cules as  long  as  the  internal  quantum  state  does  not  change.  They 
include  the  important  distortion  effect  of  the  intermolecular  potential 
on  the  incoming  and  outgoing  trajectories  which  has  been  ignored  in 
most  previous  semi-classical  calculations  and  which  is  ignored  in 
the  Born  approximation  of  quantum  mechanical  scattering  theory. 

Since  the  relative  position  of  the  colliding  molecules  is  given 
precisely  as  a function  of  time  by  these  elastic  trajectories,  the 
method  of  time  dependent  perturbations  may  be  used  to  calculate  the 
probability  that  a quantum  mechanical  transition  of  internal  coordin- 
ates occursduring  an  encounter.  When  a transition  does  occur,  it  is 
assumed  to  take  place  instantaneously,  so  that  the  direction  and  vel- 
ocity of  the  molecular  motions  change  abruptly  in  accordance  with 
♦he  classical  conservation  laws  of  energy  and  angular  momentum. 

The  trajectories  after  such  a transition  are  segments  of  other  elas- 
tic trajectories  characteristic  of  the  intermolecular  potential  deter- 
mined by  the  new  internal  quantum  states,  and  with  boundary  condi- 
tions determined  from  the  conservation  laws.  A given  internal 
transition  can  thus  lead  to  a whole  family  of  inelastic  trajectories, 
depending  where  along  the  trajectory  the  quantum  transition  takes 
place.  An  approximate  expression,  suggestion  by  the  time  dependent 
perturbation  formulae,  is  given  which  describes  the  relative  proba- 
bility of  transition  along  the  trajectory.  This  expression  is  modified 
in  Special  cases  by  a set  of  "selection  rules"  resulting  from  the  fact 
that  for  a given  transition  there  are  points  along  the  trajectory  at 
which  the  conservation  laws  cannot  be  satisfied. 

The  semi-  classical  formulation  of  the  collision  problem  has 
Esveral  advantages  over  the  quantum  mechanical  scattering  formula- 
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tion:  It  is  not  necessary  to  assume  that  the  energy  exchange  in  an 
inelastic  collision  is  small,  so  that,  for  example,  the  method  would 
be  expected  to  be  applicable  to  a study  of  rotational  energy  transfer 
in  collisions  between  hydrogen  molecules.  The  description  of  the 
relative  translational  motion  of  the  molecules,  which  leads  to  the 
greatest  difficulties  in  the  quantum  mechanical  approach  through  an 
infinite  set  of  coupled  radial  equations,  can  in  principle  always  be 
given  in  the  semi-classical  formulation  by  direct  integration  of  the 
Newtonian  equations  of  motion.  Finally,  the  semi-classical  formula- 
tion naturally  falls  into  several  distinct  steps  each  with  clear  physical 
significance,  so  that  the  significance  and  effects  of  various  simplify- 
ing approximations  are  perhaps  more  readily  seen  than  in  the  quantum 
mechanical  method. 

Following  a general  presentation  of  the  semi-classical  theory 
of  inelastic  collisions,  we  give  a detailed  description  of  the  theory 
for  the  case  of  the  transfer  of  rotational  energy  in  collisibns  between 
rigid,  linear  molecules  containing  ideal  electric  dipoles.  This  model 
is  of  particular  interest  in  kinetic  theory  since  the  high  degree  of 
asymmetry  in  the  inter  molecular  potential  would  seem  to  make  appli- 
cation of  the  Enskog  theory  for  spherical  molecules  inappropriate, 
even  as  a first  approximation.  We  include  sample  calculations  of 
the  effective  intermolecular  potential  and  probability  of  transition 
along  the  trajectory  for  typical  cases.  These  sample  calculations 
are  of  value  for  indicating  simplifications  which  might  be  made  ir.  an 
approximate  evaluation  of  the  differential  scattering  cross  sections. 
The  details  of  a sample  inelastic  trajectory  indicate  the  importance 
of  the  inelastic  transfer  on  the  trajectory. 

The  effective  intermolecular  potential  arising  from  the  adiaba- 
tic interaction  of  such  rotating  linear  dipoles  has  been  calculated  for 
several  typical  cases  from  well  known  formulae  due  to  F.  London. 
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When  the  principle  rotational  quantum  numbers  of  the  colliding  mole- 
cules differ  by  one  - the  case  of  resonance  interaction  - long  range 
attractive  and  repulsive  potentials  arise.  With  the  help  of  the  semi- 
classical  theory  it  is  shown  why  the  Enskog  theory  actually  is  able 
to  predict  the  coefficient  of  viscosity  of  polar  molecules  when  a poten- 
tial energy  function  is  used  which  includes  the  long  range  interactions. 

The  dependence  of  the  rotational  transition  probability  on  rela- 
tive velocity,  impact  parameter,  and  magnitude  of  the  energy  transfer 
is  calculated  for  some  sample  cases.  Since  the  transition  probability 
varies  rapidly  as  these  variables  are  changed,  it  is  seen  that  the 
important  contributions  to  the  cross  sections  would  be  expected  to 
come  from  collisions  characterized  by  rather  narrow  ranges  of  these 
variables. 

Calculations  of  the  relative  probability  that  transitions  occur 
at  various  points  along  a trajectory  indicate  that  the  distance  of  closest 
approach  of  the  molecules,  frequently  taken  as  the  most  probable 
place  for  a transition,  may  actually  be  a very  unlikely  place.  Instead, 
the  regions  where  the  effective  inter  molecular  potential  is  changing 
most  rapidly  with  time  are  indicated  as  the  most  likely  points  for  a 
transition.  This  fact,  too,  would  be  of  great  value  for  approximate 
calculations  of  the  cross  sections. 

The  special  case  of  the  transition  probability  when  the  mole- 
cules undergo  a resonant  transition  is  discussed.  As  long  as  the 
molecules  interact  strongly  enough  to  give  an  angle  of  deflection  appre- 
ciably different  from  zero,  it  is  shown  that  the  classical  analogue  of 
the  transition  probability  corresponds  to  a frequent  exchange  of  the 
quantum  of  rotational  energy  during  the  collision.  Thus  v.^hen  the 
molecules  finally  separate  after  the  encounter  there  is  equal  probabil- 
ity that  either  one  of  them  will  have  the  extra  quantum  of  rotational 


energy. 
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This  detailed  consideration  of  rotational  transitions  in  polar 
molecules  is  followed  by  a brief  discussion  of  the  form  the  theory 
would  take  in  the  consideration  of  vibrational  and  electronic  transitions. 


1.2  Inelastic  Collisions  and  the  Kinetic  Theory  of  Gases 

The  kinetic  theory  of  gases  composed  of  particles  which  inter- 
act according  to  a spherically  symmetric  interaction  potential  is 

completely  understood^.  That  is,  the  transport  coefficients  of  a 
4 

monatomic  gas  may  be  calculated  to  any  degree  of  accuracy  for 

an  arbitrary  iuiermolecular  potential.  The  derivation  begins  with 

5 

the  Boltzmann  equation  which  specifies  the  singlet  distribution  of 
molecules  in  phase  space.  Solutions  of  this  equation  due  to  Enskog 


4.  Actually  the  theory  may  be  used  satisfactorily  to  predict  the 

transport  coefficients  except  thennal  conductivity,  for 
molecules  like  0^  and  which  are  nearly  spherical. 
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MTGL  p.  444. 
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and  Grad  lead,  finally,  to  the  formulae  for  calculation  of  the  trans- 
port coefficients.  The  dynamical  information  required  for  these 
calculations  is  simply  the  angle  of  deflection  in  a bimolecular  colli- 
sion X as  a function  of  the  collision  parameter  b and  the 
relative  velocity  g.  (See  Fig.  2.  1). 


Fig.  2.  1 A Typical  Bimolecular  Collision  Pictured 
in  a Coordinate  System  Translating  with  the  Center  of 
\ias8  of  the  System. 


The  classical  Boltzmann  equation  does  not  apply  if  the  molecules 
have  internal  degrees  of  freedom.  By  assuming  that  the  internal 
degrees  of  freedom  could  be  treated  quantum  mechanically  while  the 
translational  motion  was  treated  with  classical  mechanics.  Wang 
Chang  and  Uhlenbeck  were  able  to  modify  the  theory  so  that  it  would 
describe  this  more  complicated  situation.  In  principle,  their  method 
amounts  to  considering  the  gas  as  a mixture,  with  each  set  of  internal 
quantum  states  defining  a separate  species.  A separate  Boltzmann 
equation  gives  the  distribution  function  for  each  species,  and  the 
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equations  are  coupled  in  the  sense  that  collisions  tranrfer  molecules 
from  one  species  to  the  other.  A solution  of  these  coupled  equations, 
similar  to  Enskog's,  then  leads  to  the  expressions  for  the  transport 
coefficients.  Wang  Chang  and  Uhlenbeck  give  the  formulae  for  the 
transport  coefficients  for  the  two  limiting  cases  of  very  easy  and 
very  difficult  exchange  of  internal  energy.  Whereas  the  theory  for 
spherical  molecules  depended  only  on  the  angle  of  deflection 
the  theory  for  polyatomic  molecules  depends  on  the  differential 

scattering  cross  sections,  The  differential 

— 

cross  section  for  a collision  process  is  defined 

as  the  fraction  of  the  molecules  with  relative  velocity  g which 
are  scattered  from  a uniform  incident  beam  of  flux  one  molecule 
per  unit  area  per  unit  time,  into  a unit  solid  angle  in  the  direction 
Xj  while  the  internal  quantum  numbers  specifying  the 

system  go  from  i and  j to  k and  1. 

Although  we  do  not  intend  to  describe  either  of  these  theories 
in  any  detail,  it  will  be  instructive  to  compare  typical  results  from 
the  two  theories  in  order  to  see  the  effect  of  the  internal  degrees  of 
freedom.  For  molecules  with  spherical  potentials,  the  first  approx- 
imation to  the  reciprocal  of  the  viscosity  coefficient  is  given  by^ 

CD  0> 


lb  V 

S iU  cw  JiT 


e h dh  (1.2- 1) 
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Here  [^j,  is  the  first  approximation  to  the  coefficient  of 
viscosity,  /M.  is  the  reduced  mass,  ^ the  relative 

speed  of  the  molecules,  b the  collision  parameter,  and 


is  the  reduced  mass, 


the  relative 


speed  of  the  molecules. 


is  the  angle  of  deflection  of  the  collision.  These  quantities  are  best 
understood  by  reference  to  Fi  g 2.  1 where  a typical  collision  is 
pictured  in  a coordinate  system  translating  with  the  center  of  mass. 

The  important  thing  to  notice  about  Eq.  1.2-1  is  that  it 
amounts,  essentially,  to  a calculation  of  the  quantity  ~ 

averaged  over  all  collisions.  Since  7L~n  h •i-b  is  the  probability 

that  a collision  in  a dilute  gas  near  equilibrium  will  occur  with 
collision  parameter  between  b and  b + db  , and  since  ^ 6.  ®'*‘ 
is  proportional  to  the  probability  that  the  relative  velocity  of  the 
encounter  will  lie  between  g and  g + dg,  the  averaging  corres- 
ponds to  an  equilibrium  distribution  of  all  possible  collisions.  The 
angle  of  deflection,  X(^  depends  only  on  b and  g , so 

the  average  may  be  accomplished  by  a double  integration  over  these 
two  initial  conditions. 

When  the  molecules  possess  internal  degrees  of  freedom,  and 

when  the  energy  of  those  degrees  of  freedom  exchanges  readily  with 

7 

the  tran'.lational  energy,  Wang  Chang  and  -hlenbeclr  show  that 
the  expji^noion  for  the  reciprocal  of  c oc  fficient  of  viscosity  is 

^ y ^-TT— Ui 

[vj  I ^ 

"'i  ' i ■>  (1  - 

O ^ _ (i. 


-M-  % 


7.  Loc.  ci!;.  p.  20 


10. 


1 


Here  to,  is  the  internal  energy  of  mclectAle  a when  it  is  in 
the  i'th  quantum  state,  is  the  net  change  in  internal  energy, 

and  Ip  is  the  azimuthal  scattering  angle  which  tells  how  far  the 
final  trajectory  is  bent  out  of  the  plane  of  the  initial  trajectories.  The 
significance  of  PC  is  unchanged. 

Equation  1.2-2  differs  from  Eq.  1.2-1  in  two  respects.  The 
property  of  the  collision  which  is  to  be  averaged  is  now 

I X f ^ ^ ^ t-t,  ( I " (1. 2-3) 


Clearly  this  reduced  to  the  expression  in  Eq.  1.2-1  when  the  collision 
is  elastic,  since  then  is  zero.  In  the  second  place,  we  now  label 

the  encounters  not  with  the  initial  conditions  b and  g,  but  with  the 
initial  velocity  g and  the  final  scattering  angles  X and  'P  f 
together  with  the  initial  and  final  internal  quantum  states.  This 
change  in  point  of  view  complicates  the  averaging  process  and  necessi- 
tates introduction  of  the  differential  scattering  cross  sections 
X X quantity  ^ ^ 

now  plays  a rule  analogous  to  , enabling  us  to  make  an  aver- 

age over  all  possible  encounters  even  though  we  are  counting  them  with 
respect  to  labels  characteristic  of  the  final  state*  . Evaluation  of 
these  cross  sections  represents  the  chief  problem  in  the  application 
of  the  kinetic  theory  of  polyatomic  molecules.  In  the  next  section  we 
describe  the  quantum  mechanical  significance  and  evaluation  of 
the  cros.T  sections. 


* This  change  of  variables  could  be  carried  out  in  the  classical 
formula  by  changing  from  an  integration  over  b to  one  over 
X then  3.77  [_;j^  . 

The  quantity  in  brackets  is  just  what  is  usually  defined  as  the 
classical  differential  scattering  cross  section.  (See  for 
example,  H.  Goldstein,  Classical  Mechanics,  Addison  Wessley, 
1951).  For  the  case  of  rigid  spheres  it  is  simply  the  geometric 
cross  section  'Y  where  d is  the  radius  of  the  sphere. 


i 

I 

i 
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I 
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1. 3 The  Quantum  Theory  of  Scattering  ’ 

A detailed  discussion  of  the  quantum  theory  of  scattering  would 
be  out  of  place  here,  but  we  shall  outline  eno\igh  of  the  theory  to 
show  the  quantum  mechanical  significance  of  the  cross  section,  and 
to  indicate  the  chief  approximate  methods  which  have  been  applied 
to  studies  of  molecular  collisions.  The  approximate  methods  will 
be  of  use  for  purposes  of  comparison  with  their  counterpartp  in 
the  semi-classical  theory  to  be  presented  in  Chapters  II  and  HI. 

Suppose  we  consider  isolated  molecules  a and  b,  vhose 
internal  Hamiltonians  are  Ho-  and  Hb  and  whose  wave 
functions  Axt,  and  M-i,  and  energy  levels  Eo.  and  Ei 


8.  Most  of  this  section  is  based  on  material  in  W F.  Mott 
and  H.  S.  W.  Massey,  The  Theory  of  Atomic  Collisions 
Oxford  Press,  seconded.  (1949). 

9.  For  recent  important  developments  in  the  fundamental 
theory  of  scattering, see  : 

D.  R.  Bates,  A.  Fundaminsky,  H.  S.  W.  Massey,  and 
J.  W.  Leech,  Phil.  Trans.  Roy.  Soc.  A-860,  93,  (1950). 

Breit,  Rev.  Mod.  Phys.  23,  238  (1951) 


S.  Altschuler  Pnys.  Rev.  1278  (1953);  95.  546  (1954) 

C.  F.  Curtiss,  J.O.  Hirschfelder  and  F.  T.  Adler 
J.  Chem.  Phys.  1^,  1638  (1950) 

C.  F.  Curtiss,  J.  Chem.  Phys.  1199  (1953) 
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are  obtained  by  solutions  of  the  equations 

LHbiw-  EJJ  O 


(1.344 


Here  I»-  and  denote  all  of  the  coordinates  necessary  to 

specify  the  configurations  of  molecules  a and  b relative  to 
their  mass  centers. 

If  position  of  molecule  b relative  to  molecule 

a,  g the  relative  velocity,  the  reduced  mass,  and 

8ai>,  ]>,£>)  the  potential  energy  of  interaction,  then  the  wave 


equation  of  the  colliding  system  in  relative  coordinates  is 


. 10 


■ Fa.'  ^ 


(i.3-2) 


r:  O 


The  problem  is  to  solve  Eq.  1.  3-2  with  boundary  conditions  cor- 
responding to  a molecular  collision.  L>et  us  look  for  the  solution 
which  describes  the  scattering  of  a uniform  beam  of  molecules  in 
state  i from  a molecule  in  state  j.  TaJcing  the  incident  beam 
along  the  positive  z axis,  and  neglecting  symmetry  effects,  the 

solution  may  be  expanded  in  the  complete  set  of  functions 
A 

M-k'V--)  (including  the  continuum,  for  ionized  states). 

Thus 


$ 


t J 


r-  > «■  / 

j (l?^Lb)  yUaCt“)  ^b(^. 


(1.3-3) 


10. 


W.F.  Mott  and  H.  S.W.  Massey,  ]oc,  cit.  p.  140. 


To  satisfy  the  boundary  conditions,  we  must  have  for  large  I 
h\  IS^»!  S (1.3-4) 


I Rab  I 


It  may  be  shown  that 


(1-3-5) 


— J/yy> 


The  radial  functions  LE<^^)  are  given  by  the  set  of 

coupled  equations 


.y.  A ^ 

Mi(Cu)  Ml T(ry)  & < &) ^ 


11.  ibid.  p.  137 
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Obtaining  solutions  to  this  infinite  set  of  coupled  radizil 

equations  constitutes  the  chief  problem  of  the  quantum  mechanical 

12 

method.  As  Zener  has  pointed  out  , it  is  the  difficulty  of  handling 
the  radial  equations  in  the  quantum  mechanical  method  that  suggests 
development  of  a semi-c lassical  theory  for  investigating  low  velocity 
molecular  collisions.  The  translational  motion  of  molecules  in 
thermal  collisions  is  very  nearly  classical,  so  that  one  might  hope 
for  a considerable  simplification  when  the  problem  is  reformulated 
in  such  a way  as  to  utilize  Newtonian  mechanics  for  a description 
of  the  relative  molecular  motion. 

We  shall  now  list  three  approximate  solution.s  of  the  scattering 
equations  which  are  of  importance  in  the  study  of  molecular  collisions., 
One  can,  to  some  extent,  attach  a physical  significance  to  the  math- 
ematical approximations  made  in  each  case, 
a)  The  Born  Approximation. 

If  on  the  right  hand  side  of  Eq.  1.  3-6  we  take 

t j c L 015  • 

- e 

(1.3-7) 


r-  ^ 

F.j  ^ 


i,  ^ ^ 


where  u is  a unit  vector  in  the  direction  of  the  incident  beam 

0 

and  n is  a unit  vector  in  the  direction  of  then  it  may 

1 7^ 

be  shown  that 


12. 


C.  Zener,  Phys.  Rev.  557  (1931) 

W.  F.  Mott,  and  H.  S.  W.  Massey,  loc.  cit.  p.  143. 


13. 
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K( 


K 


J 


t -K«<  • ^O.b  

i’, 


(1.3-8) 


To  this  approximation,  the  distortion  of  the  incoming  and 
outgoing  waves  by  the  interaction  potential  is  minimized,  since 
they  are  taken  to  be  plane  waves,  in  a semi-classical  treatment 
this  would  correspond  to  taking  these  trajectories  as  straight  lines 
along  which  the  molecules  travelled  at  constant  velocities  character- 
istic of  the  initial  and  final  translational  energy.  Thus  the  method 
is  best  suited  to  a study  of  high  energy  collisions  where  the  trans- 
lational motion  is  only  slightly  perturbed  by  the  collision.  Because 

14 

of  the  relative  simplicity  of  the  method,  however,  Kerner  has 
used  it  for  a study  of  both  vibrational  and  rotational  energy  transfer 
in  thermal  molecular  collisions. 

b)  Distorted  waves. 

If  we  make  the  less  drastic  assumption  that  the  non  diagonal 
matrix  elements  on  the  right  hand  side  of  Eq.  1.  3-6  are  so  small 
that  we  may  neglect  all  terms  except  the  diagonal  ones; 

and  the  off  diagonal  ones  connected  with  the  incident  state: 


$ 


^ ti 

1“  • 

• <-  j 


(^Ra.bj 


14. 


E.  H.  Kerner,  Phys.  Rev.  92,  899  (1953  and  91  , 

1174.  (1953). 
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then  it  may  be  shown 


15 


that 


(1.3-9) 


Here 


is  the  solution  of 


[v" . K.r  - F:/m  = o 


{1.  3-10) 


and  the 


Kf 

J'ci  ( ll^») 


are  solutions  of 


[_V"  ^ K„<’ 


V _ ^ /U. 


(1.3-11) 


In  this  approximation  some  of  the  distortion  of  the  incident  and 
outgoing  waves  by  the  scattering  field  is  taken  into  account.  Instead 
of  using  plane  waves  for  the  incident  and  scattered  waves,  we  now 


15. 


W.  I*.  Mott  and  H.  S.W.  Massey,  loc.  cit.  p.  146. 
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have  waves  distorted  by  the  diagonal  elements  of  the  interaction 
potential.  This  is  the  method  most  frequently  used  for  the  investi- 
gation of  molecular  collisions  . 

c)  Perturbed  stationary  states. 

In  this  method  the  original  derivation  is  varied  slightly.  It  i^ 
assumed  that  the  encounter  is  nearly  adiabatic,  and  that  for  every 
inter  molecular  distance  the  equation 


Q*  -r,.  BsLb)  - E i giibj  *4^  I gab)  - ^ 


(1.3-12) 


has  been  solved.  The  tp  I Rat.)  form  a complete 

set  so  that  the  exact  solution  may  be  expanded  in  terms  of  these 
functions  instead  of  the  products 
Thus  the  solution  is  expanded  in  the  form 


.4  (1-3-13) 

The  derivation  is  carried  through  with  the  same  approximations 
as  for  the  method  of  distorted  waves.  If  in  the  final  result,  the 


16.  C.  Zener,  Phys.  Rev.  j[7,  556  (1931). 

17.  R.  S.  Roy  and  M.  E.  Rose,  Proc.  Roy  Soc.  A 149,  511  (1935) 

18.  R.  N.  Schwartz,  Z.  I.  Slawsky,  and  K.  F.  Hersfeld, 

J.  Chem.  Phys.  1591  (1952). 

19.  K.  Takayanagi,  Prog.  Theoret.  Phys.  111  (1952). 

20.  J.  C.  Beckerle,  J.  Chem.  Phys.  21,  2034  (1953). 

21.  R.  N.  Schwartz  and  K.  F.  Herzfeld,  J.  Chem.  Phys.  22, 

767  (1954). 

J.  C.  Broot,  J.  Chem.  Phys.  22,  1189  (1954.). 
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gAb)  are  replaced  by  their  first  order  pertur- 
bation approximation  in  terms  of  the  JJ.ou(Xs-)  and  fa 

.23 


there  is  obtained 


-rKf  h-  Jill 

K., 


A 


(1.3-14) 


From  this  approximate  expression,  it  may  be  seen  that  the 
result  is  similar  to  the  distorted  wave  approximation,  with  full 
account  taken  of  the  distortion  effect  on  the  incoming  and  outgoing 
trajectories.  Although  the  method  appears  well  suited  to  the  study 
of  molecular  collisions,  it  ia  difficult  to  apply  since  it  requires 
a reasonably  exact  knowledge  of  the  adiabatic  functions.  The  exci- 
tation of  He  atoms  by  protons,  and  the  electron  capture  from 

24 

He  by  protons  were  studied  some  years  ago  by  this  method. 

Interest  in  the  method  has  recently  been  revived  and  semi-classical 

formulae  for  the  total  transition  probability  have  been  discussed 

25 

in  terms  of  the  method  by  Bates  et  al.  It  will  become  evident 

in  Chapter  II  that  our  semi-classical  description  of  the  total  colli- 
sion process  bears  a close  resemblance  to  this  method. 


23.  W.  F.  Mott  and  H.  S.  W.  Massey,  loc.  cit.  p.  156. 

24.  Massey  and  Smith,  Proc.  Roy.  Soc.  A 142 , 142  (1933). 

25.  D.  R.  Bates,  H.  S.  W.  Massey,  and  A . L.  Stewart, 

Proc.  Roy.  Soc.  A 216,  437  (1953). 


19. 
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1.4  History  of  the  Semi-Clas sical  Method  for  Studying 
Collision  Problems.  

Shortly  after  the  introduction  of  quantum  mechanical  methods 
for  studying  inelastic  collision  phenomena,  workers  recognized  the 
possibility  of  combining  a classical  description  of  the  translational 
motion  with  a quantum  description  of  the  internal  motions.  In  all 
of  these  studies  the  aim  was  to  calculate  the  total  probability  of  some 
change  in  the  internal  state,  such  as  electronic  excitation  of  ionization. 
The  details  of  the  trajectories  in  such  studies  were  of  secondary  im- 
portance, in  that  the  classical  trajectory  was  used  only  to  give  the 
explicit  time  dependence  of  an  interaction  potential.  This  quantity 
is  not  sensitive  to  slight  variations  in  the  curvature  of  an  orbit,  so 
that  in  most  cases  straight  line  trajectories  were  found  to  be  adequate 
for  the  problems  considered.  For  application  to  the  kinetic  theory  of 
gases,  we  require  much  more  detailed  and  more  accurate  information 
about  the  shape  of  the  trajectories.  Since  oxir  theory  is  based  on 
these  earlier  works,  we  shall  review  briefly  some  of  the  high  points 

in  the  development  of  the  semi-classical  method. 

26 

Faxen  and  Holtzmark  were  perhaps  the  first  to  emphasize 
the  fact  that  the  semi-classical  method  depended  on  a quantum  mechan- 
ical separation  of  internal  coordinates  from  the  coordinates  of  relative 
motion.  Any  use  of  the  semi-classical  method  must  begin  with  a 

27 

justification  of  this  separation  for  the  particular  problem.  Gaunt 

made  one  of  the  first  calculations  of  a collision  problem  by  the  semi- 

classical  method  when  he  considered  the  scattering  of  alpha  particles. 

Z8 

His  results  were  quite  poor,  but  it  was  later  shown  by  Mott  that 

26.  Faxen  and  Holtzmark,  Zeit.  f.  Physik  45,  311  (1927). 

27.  Gaunt,  Proc.  Camb.  Phil.  Soc . , 23,  732  (1927). 

28.  N.  F.  Mott,  Proc.  Camb.  Phil.  Soc.  27,  553  (1931). 


20. 

the  bad  results  were  due  to  the  approximation  of  some  integrals  rather 

than  any  inadequacy  of  the  nnethod.  In  the  same  paper,  Mott  put  the 

semi-classical  pr  .cedure  on  more  firm  ground  by  demonstrating  its 

relation  to  the  Born  approximation  of  quantum  mechanical  scattering 

theory,  and  it  was  he  who  called  it  the  method  of  impact  parameters. 

29 

Frame  calculated  the  probability  of  atomic  excitation  by  fast  alpha 

particles,  using  both  the  Born  approximation  and  the  semi-classical 

method.  The  agreement  was  good. 

The  first  application  of  the  method  to  inelastic  molecular  colli- 

30 

sions  was  made  by  Zener  when  he  considered  vibrational  excitation 

produced  by  head  on  molecular  collisions,  and  showed  that  the  method 

was  considerably  more  simple  than  the  equivalent  treatment  by  the 

31 

method  of  distorted  waves.  Recently  Widom  and  Bauer  have  rcviv'»d 
Interest  in  the  method  by  extending  Zener's  treatment  to  include 

25 

collisions  in  which  the  im.pact  parameter  differs  from  zero.  Bates 
and  coworkers  have  stressed  the  importance  of  the  distortion  of 
trajectories  in  slow  collisions  and  have  extended  the  method  in  this 
direction  in  connection  with  their  work  on  the  quantum  mechanical 
method  of  perturbed  stationary  states. 


29.  J.  W.  Frame.  Proc.  Camb.  Phil.  Soc.  Z7,  511(1931) 

30.  C.  Zener,  Phys.  Rev.  38,  277  (1931) 

B.  Widom  and  S.  H.  Bauer,  J.  Chem.  Phys.  1670  (1953) 


31. 
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U.  THE  SEMI-CLASSICAL  DESCRIPTION 
OF  A MOLECULAR  COLLISION 


In  this  chapter  we  present  the  semi-classical  description  of  a 
general  molecular  collision.  In  section  2.  1 we  review  the  starting 
equations  and  present  the  basic  assumptions  of  the  method.  Next  we 
apply  the  equations  to  a description  of  collisions  in  which  no  internal 
transition  takes  place.  In  section  2.  2 we  extend  the  method  to  include 
a description  of  the  internal  wave  functions  and  trajectories  when  a 
transition  of  internal  coordinates  takes  place  during  the  encounter. 


2.1  Starting  Equations  and  Elastic  Encounters 

Consider  a system  of  two  colliding  molecules,  "a"  and  "b", 
with  masses  m^  and  m^.  The  motion  of  the  system  may  be  des- 
cribed in  terms  of  that  of  the  center  of  mass  of  the  complete  system, 
the  relative  motion  of  the  centers  of  mass  of  the  two  bodies  ""^d  the 
motion  of  the  individual  molecules  relative  to  their  own  centers  of 
mass.  Of  there,  the  motion  of  the  center  of  mass  of  the  entire  system 
is  irrelevruit  and  can  always  be  separated  out  and  the  Schreedinger 
equation  describing  the  two  molecule  system  may  be  written  in 
form  . 


1 


N.  F.  Mott  and  II.  S.  W.  Massey,  loc . cit.  p.  139. 
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-|i  " Er*J@(fe,n..  Sx^;  = O (2.  i-i) 


Here 


and  r , denote  all  of  the  coordinates  necessary  to 
— b 


specify  the  configurations  of  molecules  a and  b relative  to 

their  own  mass  centers,  and  H and  H,  are  the  Hamiltonians 

a b 

of  molecules  a and  b when  they  are  isolated  in  space.  R , 

— ab 

gives  the  orientation  of  the  center  of  mass  of  molecule  b relative 
to  that  of  molecule  a,  viewed  in  a coordinate  system  translating 
with  the  center  of  mass  of  the  whole  system,  and  the  total 

energy  of  the  system  relative  to  the  center  of  mass.  The  quantity 
(&  G>  potential  energy  of  interaction  of 

the  molecules,  and  the  reduced  mass  of  the  system 
is  given  by 

/We.  /Wfe 

h-  = 


/rria,  + /Wb 

The  internal  Hamiltonians  H 


and  H^  may  be  used  to  define 


a set  of  eigenfunctions  XjUo.  , M,t  and  energies  ta,  , tZt> 
describing  the  motions  of  the  isolated  molecules 

- El]  - O 


(2.  1-2) 


The  fact  that  only  the  relative  separation  of  the  two  molecules 

appears  in  Eq.  2.2-1  indi  cates  that  the  relative  motion  of  the  mole- 

2 

cules  may  be  viewed  in  terms  of  an  equivalent  one  particle  problem  . 


2 


MTGL.  p.  49- 
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The  information  describing  the  motion  of  the  two  molecules  in  the  cen- 
ter of  mass  coordinate  system  may  be  obtained  by  imagining  that  one 
molecule  is  fixed  at  the  origin  and  that  the  other  movea  relative  to  it 
with  the  same  inter  molecular  potential  as  before,  providing  only  that 
the  moving  molecule  is  assumed  to  have  a mass  equal  to  the  reduced 
mass  of  the  two  molecule  system.  The  description  of  the  internal  mo- 
tions of  the  two  molecules  is  not  altered  by  this  changed  point  of  view, 
since  the  relative  coordinate  system  translates  uniformly  with  the 
center  of  mass,  but  does  not  rotate  with  respect  to  axes  fixed  in  space. 
Because  of  the  conceptual  simplicity  of  the  one  pa'^ticle  description  of 
the  encounter,  we  shall  use  it  extensively  in  the  rest  of  the  discussion. 

For  the  semi-classical  solution  of  Eq.  2.  1-1  we  assume  that 

3 4 

a Born  - Oppenheimer  type  separation  ’ of  the  internal  coordinates 
from  the  coordinates  of  relative  motion  may  be  used.  That  is,  the 
solutions  of  Eq.  2.  1-1  are  assumed  to  be  expressable  as  a product 

O ([<!-,  r»  / f ti-  I & I E'-**)  X ( 

where  and  are  given  by  the  equations 


+nb|L'»)  - rt.,Rab)  (2.1-3) 

and 


A (dav 


(2.1  -4) 


The  wave  function  not  to  be  confused  with  the  azimuthal 

angle  of  scattering 

Equation  2.  1-3  gives  the  wave  functions  and  energies  of  the 
hypothetical  situation  in  which  the  relative  orientation  of  the 


3.  M.  Born  and  J.  R.  Oppenheimer,  Ann.  Phys.  84,  457  (1927) 

4.  Derivations  of  the  separation  and  ways  of  determining  the 
error  are  further  discussed  in  H.  Pelzer  and  E.  P.  Wigner. 
Z.  Physik.  Chem.  B 15  , 445  (1932)  and  MTG.L  p.  925. 
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molecules  is  held  fixed.  The  wave  functions  ^ 

and  energies  of  this  solution  thus  depend  only  parame- 
trically on  the  intermolecular  orientation  Ra,b  The  quantum 

numbers  i and  j represent  all  necessary  quantum  numbers  for 


specifying  the  internal  state.  A representation  is  chosen  such  that 
functions  as 


becomes  a simple  product  of  the  separate  molecule  wave 


R is  increased  to  infinity.  Thus 

^ ‘ ® ^ (C?J 

lii  ^ 

YirJt 


(2.  1-5) 


For  this  reason  we  refer  to  the  y "adiabatic 

clamped"  solutions  of  the  internal  motions.  The  energies 
define  a set  of  potential  energy  surfaces,  one  for  each  set  of  quantum 
numbers,  which  govern  the  relative  motion  of  the  moecules.  Equation 
2.  1-4  then  gives  the  quantum  mechanical  description  of  the  relative 
motion  of  the  molecules  on  such  a potential  energy  surface.  It  is 
convenient  to  define  two  new  quantities  and 


such  that 


( J 


This  is  reasonable  and  useful  since  now 


(2.  1-6) 


I 


vanishes 


CP  (Rat) 

as  j Rg.b|  goes  to  infinity  and  is  indeed  the  effective  intermolecular 

» r~  ^ ^ • 

potention,  and  relative  translational  kinetic  energy 

of  the  colliding  system  when  the  encounter  begins. 


E- 


=-  ± 


With  this  notation  Eq.  2.  1-4  becomes 


H. 
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The  semi-classical  formulation  of  the  problem  is  completed  by 

assuming  that  £q.  ?.  1-7  for  the  relative  motion  of  the  molecules 

may  be  solved  in  the  classical  limit. 

Now  we  are  in  a position  to  give  the  semi-classical  description 

of  the  elastic  collision  between  molecules  which  start  out  in  the 

states  u^  and  u^  . The  trajectory  is  given  by  the  motion  of 

° ^ <:d 

a particle  of  mass  moving  in  the  potential  field  ^ 

with  initial  conditions  of  energy  and  angular  momentum  given  by  the 

initial  relative  velocity  g and  collision  parameter  b . The 

wave  functions  describing  the  internal  state  of  the  molecules  is  given 

by  the  adiabatic  clamped  functions  evaluated  along 

the  trajectory. 


2.  2 Inelastic  Encotmters*' 

An  inelastic  encoimter  is  described  by  assuming  that  a transition 
of  the  internal  quantum  state,  say 

D>i  RfO 

takes  place  at  some  point  along  the  trajectory.  The  subsequent  rela- 
tive motion  is  described  by  motion  on  the  new  potential  energy  surface 
with  boundary  conditions  determined  by  application  of 
the  classical  conservation  laws  of  energy  and  angular  momentum  at 


5.  Kallman  and  London,  Z.  Phys.  Chem.  B 2 , 207  (1929)  give  a 

thorough  discussion  of  the  many  types  of  non  adiabatic  molecular 
collis  ions . 
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the  point  of  transition.  In  piirt  ^ of  this  section  the  probability  of 
transition  is  discussed  in  terms  of  time  dependent  perturbation  theory, 
and  an  approximate  expression  is  presented  for  determining  the  rela- 
tive probability  that  the  transition  took  place  at  various  points  along 
the  trajectory.  In  part  the  use  of  the  classical  conservation  laws 

for  determining  the  new  trajectory  is  discussed. 


a)  The  internal  transition. 

The  classical  elastic  trajectories  give  the  inter  molecular 

separation  R as  an  explicit  function  of  time,  so  that  the  method 

of  time  dependent  perturbations  is  indicated  for  an  investigation  of 

the  probability  of  transition.  One  slight  modification  of  the  theory  is 

necessary  for  our  purposes,  however.  Ordinarily  one  takes  for 

unperturbed  wave  functions  the  simple  product  functions  corresponding 

6 

to  the  non  interacting  system.  Then  one  finds  that  the  probability 
that  a system  which  at  time  t = ~CD  was  in  the  pure  initial  state 
whose  qxiantum  nunnbers  are  specified  by  i and  j is  found  in  the 
final  state  specified  by  k and  1 at  time  t , is  given  by 


' ij 


TT  e 


-00 


(2.2-1) 


j j I 

Here  E , E , El  , and  E,  are  the  internal  energies  of  the 
a a b b 

isolated  molecules  as  defined  in  Eq.  2.  1-2,  and  the  matrix  element 
in  the  inter  grand  is  given  by 


6.  P,  A.  M.  Dirac,  Froc.  Roy.  Soc.  A112  . 661  (1926)  and 
A 114  . 243,  (1927) 

A derivation  of  the  method  is  included  in  most  elementary 
quantum  mechanics  texts,  see  for  example,  L.  I.  Schiff, 
Quantum  Mechanics,  McGraw  Hill,  1949,  p.  189. 
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We  suppose  that  "t  -O  is  taken  at  the  distance  of  closest  approach 
of  the  two  molecules. 

In  our  formulation  of  the  collision  problem  in  Sec.  2.  1 we 
assumed  that  the  encounters  were  nearly  adiabatic,  so  that  the  adia- 
batic functions  obtained  by  solution  of  Eq.  2.  1-3 

provide  a more  accurate  description  of  the  internal  coordinates  during 
the  encounter  than  do  the  • Taking  these  adiabatic  functions 

(which  depend  on  time  through  ^ continuously  changing  set 

of  unperturbed  functions,  it  is  shown  in  appendix  I that  the  transition 
probability  ia  given  by 


P''w 


I Rok(t';  - 


(2.  2-2) 


Here  the  superscript  (1)  on  the  transition  probability  indicates 
that  the  adiabatic  wave  functions  and  energies  are  taken  for  the  unper- 
turbed set.  Note  that  now  the  time  derivative  of  the  interaction  poten- 
tial appears  explicitly,  and  that  the  energies  of  the  unperturbed  states 

vary  with  the  time  through  R , . 

— ab 

The  most  interesting  feature  of  this  form  of  the  transition 
probability  is  the  appearance  of  the  time  rate  of  change  of  the 
potential,  rather  than  the  potential  itself,  as  in  the  expression  for 


In  the  limit  of  an  infinitely  dow  encounter  we  expect 


no  transitions  at  all,  for  although  the  wave  functions  are  distorted  by 
the  interaction  potential  during  such  an  encounter,  they  return  adia- 
batically  to  their  initial  state  as  the  interaction  decreases. 
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That  this  is  indeed  the  case  may  be  seen  by  writing  Eq.  2.  2-2 
as  a function  of  the  arc  length  along  the  traject^^ry;  s,  instead  of  as 
a function  of  the  time.  Then 


tij 

‘ ! »Kt 
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(2.  2-2b) 


Now,  as  the  relative  velocity  vanishes,  a small  change  in  s corres- 
ponds to  a large  change  in  the  time  t(s)  so  that  the  exponential  term 
oscillates  more  and  more  rapidly  as  the  relative  velocity  is  decreased. 
Thus  successive  contributions  of  the  integrand  tend  to  cancel,  and  it 
may  be  seen  at  least  qualitatively  that  as  the  relative  velocity  vanishes, 
the  integral  goes  to  zero. 

Thus  the  adiabatic  distortion  of  the  wave  functions  has  been 
separated  out,  and  v.'e  are  looking  only  for  transitions  among  such 
distorted  wave  functions,  which  are  caused  by  the  relative  velocity. 
This  may  be  visualized  as  transitions  among  the  distorted  potential 
energy  surfaces  of  the  system,  as  in  Fig.  2.2-1. 
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Fig.  2.2-1  A quantum  transition  of  internal  coordinates 
viewed  as  a jvimp  between  potential  energy  surfaces. 
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Here  we  picture  schematically  several  potential  energy  surfaces  of 

the  system  during  an  encounter  as  a function  of  time.  The  function 

I' gives  the  probability  that  a system  which  started  out 

on  the  ij  energy  surface  will  at  time  t have  jumped  across  to 

7 

the  kl  energy  surface.  Bohm  points  out  that  this  formulation 
of  the  time  dependent  perturbation  in  terms  of  adiabatic  fu  nctions 
allows  one  to  consider  systems  in  which  the  perturbation  potential  is 
large,  as  long  as  the  rate  of  change  of  the  potential  is  small. 

In  most  collision  problems,  one  asks  about  the  total  probability 
that  a transition  has  taken  place  due  to  a perturbation  during  the 
whole  course  of  the  collision;  that  is,  one  is  interested  in  the  value 

III 

of  We,  however,  shall  be  interested  in  the  addi- 

tional question  as  to  the  relative  probability  that  the  transition  took 
place  in  given  short  time  intervals  during  the  time  the  perturbation 
was  acting.  From  Eq.  2.  .2-2  we  might  say  that  the  probability 
that  the  transition  took  place  between  time  t and  t + <Tt  is 
just 

(<i 

A P 
n.  1 

Jit 

but  in  quantum  mechanics  we  are  limited  in  the  kind  of  questions  we 

g 

may  ask  by  the  uncertainty  principle.  In  one  form  this  says  that 
the  product  of  the  uncertainty  in  a measurement  of  the  energy  of  a 
system,  5” E,  and  the  uncertainty  in  the  time  at  which  the  measure- 
ment was  made,  5 1 , are  related  by  the  expression 

7.  D.  Bohm,  Quantum  Theory,  Prentice -Hall,  1951,  p.  496. 

8.  L.  I.  Schiff,  loc.  cit.  p.  7. 
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In  our  transition  probability  we  are  asking  about  two  states  of  the 
system  differing  in  energy  by  an  amount 

It  t _ » 1 

5E  - 

hence  we  cannot  expect  the  expression 


cL  t 


to  give  meaningful  results  for  tinoes  less  than  about 


t 


C|^  C 

In  Sec,  3.2  where  we  calculate  sample  curves  of  and 

special  case  of  rigid,  linear  molecules  with 
ideal  electric  dipoles,  we  find  that  this  is  indeed  the  case,  and  that 
these  curves  oscillate  with  a period  equal  to  while  going 


from  zero  at  i--C0  to  the  final  value  Pti»> 

dK';js, 


at  t - CO . 


Such 


oscillations 


ke 


d.t 


negative  at  some  points  during  the 


encounter,  which  interpreted  naively  would  indicate  a negative  proba- 
bility of  transition  in  those  regions. 

Since  we  are  allowed,  quantum  mechanically  speaking,  to  ask 
about  the  probability  that  a transition  has  taken  place  in  a time  interval 
corresponding  to  one  such  complete  oscillation,  one  might  expect  to 
construct  a classical  analogue  of  the  transition  probability  by  drawing 
a smooth  curve  through  points  on  the  oscillating  curve  separated  by 
times 

I 


- 


- E 


This  was  attempted,  but  for  the  sample  curves  calculated^  the 
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procedure  was  not  satisfactory  in  that  the  resulting  curve  for  the 
integrated  transition  probability  still  oscillated  and  passed  through  a 
maximum  in  the  vicinity  of  the  distance  of  closest  approach.  ( It  is 
difficult  to  interpret  a decrease  in  the  integrated  probability  after  it 
has  passed  through  a maximum)  In  any  case  we  were  unable  to 
obtain  satisfactory  criteria  for  constructing  a classically  satisfactory 
approximation  lo  the  curve  i 1^!  by  considerations  of  this 


As  a next  approach,  one  might  try  to  obtain  a series  expansion 
of  the  transition  probability  or  of  its  derivative.  The  transition 
probability  is  approximately  of  the  form 


Pc-a)  - 


(2.2-3) 


Ehfferentiating  with  respect  to  y , this  may  be  written 


(2.2-4) 


Integrating  by  parts  repeatedly  and  noting  that  all  derivatives  of  the 
function  f{^  ) vanish  as  ^ we  obtain 


d % 
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(2.2-5) 


since  f(  ^ ) for  a dipole-dipole  interaction  is  of  the  form 


each  of  the  products  of  this  series  is  an  odd  function  of  ^ so 

that  approximating  by  any  finite  number  of  terms  in  the 

series  gives  the  unsatisfactory  result  that 
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In  the  absence  of  a satisfactory  derivation  of  the  classical  ana- 


logue  of  ra-^Kt 


serve  as  a satisfactory  appruximation 


we  suggest  that  the  following  function  may 
. ..  9 


I . C.C;  r. 
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(2.  2-6) 


This  function  has  the  following  desirable  properties; 

^01 

1.  It  is  zero  at  t = -00  and  equal  to  at  t = i- OP 

in  agreement  with  the  quantum  mechanical  formula. 

2.  It  is  monotonic,  increasing  over  the  entire  range,  corresponding 
to  the  classical  picture  that  if  transitions  are  rare,  there  is  negligible 
chance  of  multiple  or  reverse  transitions,  so  that  there  ought  to  be  a 
steadily  increasing  probability  of  finding  the  system  in  the  final 

state . 

3.  If  the  tlir4e  rate  of  change  of  the  interaction  potential  is  either 
symmetric  or  antisymmetric  about  t = 0,  Eq.  2.2-6  gives 


p it.o)  = P P (6 


Since  elastic  trajectories  are  symmetric  in  shape  about  the  distance 
of  closest  approach, 

d d 


The  classical  analogue  of  the  transition  probability  might  profita- 
bly be  investigated  from  the  point  of  view  of  Wlgner's  distribution 
function  formulation  of  quantum  statistical  mechanics,  but  no 
work  has  been  done  in  this  direction.  This  method  is  described 
in  E . P . W i gner,  Phvs.  Rev.  479  (1932)  and  in  J.  H.  Irving 
and  R.  W.  Zwanzig,  J.  Chera.  Phys.  19,  1173,  (1951). 
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and  Eq.  2.  2-6  says  that  there  is  equal  probability  of  the  transition 
occurring  on  the  approaching  and  the  receding  parts  of  a trajectory. 


This  is  justified  as  a first  approximation  by  noting  that  due  to  the 


is  unchanged  if 


P 

absolute  square  operation, 

d 

is  replaced  by  ~ • That  is,  to  a first  approximation,  the 

quantum  mechatiics  does  not  care  about  the  actual  sign  of  this  quan- 


tity, and  we  are  led  to  expect  equivalent  behavior  on  the  two  halves 
of  a symmetric  trajectory 


4.  The  transition  probability  increases  most  rapidly  in  regions 
where  the  perturbation  potential  is  changing  most  rapidly  in  time. 


This  is  consistent  with  the  adiabatic  theorem,  and  the  quadratic 

d- 

dependence  on  is  suggested  by  the  absolute  square 


5.  The  factor  ^ 


cLR 


operation  in  r.%)  . 

in  the  integrand  makes  f (t) 

(jt  i.)  ^ 

agree  with  the  general  form  of  when  they  are  considered 


as  functions  of  the  intermolecular  distance  R rather  than  the 
time^^.  That  is,  it  P(ti  is  expressed  as  an  integral 

over  R , it  takes  the  form 


10.  The  other  function  one  might  guess, 

ptr 


is  unsatisfactory  from  this  point  of  view,  since  it  gives 

11.  Omission  of  the  factor  1^  in  the  definition  of  f.'t.i,  in 

Eq.  2.  2-7,  thus  giving  qualitative  agreement  with  the  form  of 


expressed  as  an  integral  over  the  time  would  have  re- 
sulted in  an  anomcly  whe  II'.;  was  considered  as  an  integral 
over  R.  That  is,  we  would  have  obtained 


Since  at  the  distance  of  closest  approach  ^ ~ ^ this  would 
imply  a vanishing  probability  that  the  transition  took  place  be- 
tween R and  R+  dR  about  the  distance  of  closest  approcich. 
The  quantum  mechanical  formula,  Eq.  2.  2-8  , indicates  that 
this  is  actually  a very  likely  region  for  the  transition. 
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(2.  2-8) 


the  form 
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(2.2-9) 

Here  the  upper  limit  of  the  integrations  indicates  an  integration  along 
the  trajectory:  since  in  general  R passes  through  a minimum,  then 
increases  to  CC  , the  integral  will  have  two  branches  when  the  upper 
limit  corresponds  to  a point  on  the  trajectory  beyond  the  distance  of 
closest  approach. 


are  plotted  for  the  special  case  of  rigid,  linear  polar  molecules,  and 
are  compared  with  the  functions  and  , 

This  completes  the  semi-classical  description  of  the  internal 
transition  on  an  inelastic  trajectory.  The  expression 

J r 7 , , 


4lP’^]dHr 


gives  the  probability  that  the  transition  takes  place  between  time  t 
and  t + dt  during  the  encounter,  and 


cAmo 

Cpo) 


- Fs 


gives  the  total  probability  that  the  transition  takes  place  during 
the  encounter. 


b)  The  inelastic  trajectory. 

If  we  suppose  that  at  some  point  along  the  trajectory  the  internal 
wave  function  changes  from  to  , 

then  the  relative  motion  of  the  molecules  in  this  region  reflects  the 
sudden  shift  from  the  potential  energy  surface  ^ (B.^)  to  the 
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new  surface  y (RoJ^)  • The  boundary  conditions  of  energy 

and  angular  momentum  necessary  for  a precise  description  of  the 
motion  in  the  new  potential  field  Cp  are  deduced  from 

the  classical  conservation  laws  of  energy  and  angular  momentum. 

Let  us  denote  the  equivalent  one  particle  trajectory  before  the 

ii  kl 

transition  by  T *'  and  that  after  the  transition  by  T The 

kl 

quantities  necessary  for  specifying  the  initial  conditions  of  T 
may  be  defined  with  reference  to  Fig.  2.  2-2  where  we  picture  the 
details  of  the  trajectory  at  the  moment  of  inelastic  transfer. 


Rt,  J 


^ // 
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Fig.  2.  2-2  Details  of  a trajectory  at  the  moment  of  inelastic  transfer. 
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Since  we  are  required  simply  to  determine  the  orientation  of  T 
relative  to  T^'^,  the  coordinate  system  we  use  at  this  point  is  quite 
arbitrary,  as  long  as  it  is  defined  in  terms  of  the  trajectory  before 


ij 


the  transition,  T 

Suppose  the  transition  takes  place  at  point  B where  IRsbl  - Rat 


ii  kl 

Let  AB  and  BC  be  tangents  to  T and  T 


at  R 


ab’ 


and  let 


AB  lie  in  the  xy  plane,  parallel  to  the  x axis.  Let  the  perpendi- 
cular distances  of  AB  and  BC  from  the  origin  be  ^ {R‘ta)  and 


Note  that 


/C; 


is  an  instantaneous  analogue  of  the 


collision  parameter  b,  in  that  for  large  values  of  ~ b . 

Let  the  origin  and  line  BC  determine  the  plane  F which  inter 
sects  the  xy  plane  along  OB  = ^ makes  an  angle  f 

with  the  xy  plane.  Finally,  let  > 

with  the  x axis.  In  case  the  effective  potentials  a^re 

will  lie  entirely  in  the 


be  the  angle  which  R*  makes 


spherically  symmetric,  the  trajectory  T 
kl 

xy  plane  and  T will  lie  entirely  in  the  F plane.  In  this  case 
it  would  be  convenient  to  fix  the  coordinate  system  by  making  the  x 
axis  parallel  to  the  initial  asymptote  of 

kl 

The  initial  conditions  of  the  new  trajectory  T may  now  be 
specified  by  three  quantities:  giving  the  inclination  of  the 

F plane,  ^6"  locating  the  trajectory  in  the  F plane,  and 

A/ 

giving  the  new  relative  velocity.  These  three  quantities 
may  be  determined  by  application  of  the  classical  conservation  laws 
of  energy  and  angular  momentum  applied  at  the  point  of  transition. 
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1.  The  relative  velocity  after  the  transition; 

In  terms  of  quantities  already  defined,  the  conservation  of  energy 
relative  to  the  center  of  mass  of  the  whole  system  may  be  written 

Erei  - ^CS^b)  -h  ^/i[V 

(2.2-10) 

- i-  ^ ^ + (pVsai,)  iRob)] 

From  Eq  . 2.1-6 

Ea  + Eb  -t-  - Ere/Cloi^) 


so  the  first  three  terms  in  the  energy  expression  represent  the  total 
internal  energy  of  the  adiabatic  clamped  system.  The  last  term 
containing  the  relative  velocity  v^'^(R*^)  , gives  the  kinetic  energy 
relative  to  the  moving  center  of  mass  of  the  two  molec\ile  system. 
Eq.  2.2-lOmay  be  solved  immediately  for  the  relative  velocity 
after  the  transition: 


(2.  2-11) 


In  case  this  velocity  is  imaginary,  we  conclude  that  the  transition  is 
energetically  forbidden,  so  this  result  constitutes  a sort  of  "selection" 
on  the  transition. 


2.  Orientation  of  the  trajectory  in  plane  F; 

The  conservation  of  angular  momentum  may  be  expressed  as 


L ^ (isb)  ~ <4^  ■*'  ^ 


(2.  2-12) 
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Here  ^ and  are  the  internal  angular  momenta 

associated  with  molecules  a and  b before  and  after  the  transition, 
and  and  are  the  angular  momenta  associated 

with  the  relative  translational  motion  of  the  molecules.  For  example 

-f 

Equation  2.  2-12  may  be  used  to  determine  ^ L as 

follows;  By  definition, 

So  if  we  assume  for  the  moment  th».t  we  are  able  to  solve  Eq.  2.  2-12 
for  then  we  have,  simply, 

" ~7TT?5r" 

The  significance  of  this  result  as  far  as  determining  the  orientation 
of  BC  in  plane  F may  be  seen  by  reference  to  Fig.  2.2-3  where 
we  picture  the  situation  in  plane  F. 


Fig.  2.  2-3  The  trajectory  in  the  new 
plane  immediately  after  the  transition. 
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From,  the  definition  of  BC  is  to  be  constructed  so  as  to  be 

tangent  to  the  circle  of  radius  ( Rfl.b)  centered  at  the  origin 

and  lying  in  plane  F.  From  Fig.  2.  2-3  it  is  apparent  that  if  the  point 
B lies  outside  the  circle  there  are  two  such  possible  tangents,  both 
of  which  give  the  correct  translational  angular  momentum.  In  case 

/j'*-  ^ » 

i®  greater  than  ®°  that  point  B lies  inside  the 

circle,  no  satisfactory  trajectory  can  be  constructed  at  point  B 
which  will  be  consistent  with  the  conservation  of  angular  momentum. 
This  constitutes  a second  selection  rule  on  the  transition  at  point  B. 

When  the  transition  is  allowed,  the  choice  between  the  two 


possible  paths  is  really  unique.  This  may  be  argued  as  follows; 

Let  us  suppose  that  the  change  in  angular  momentum  in  the  transition 
specified  by  ij  kl  approaches  zero. 

In  the  limit  of  no  transition  a diagram  like  Fig.  2.  2-2  may  still 
be  constructed  with  the  two  possible  trajectories,  but  now  the  equa- 
tions of  motion  clearly  require  that  the  particle  continue  on  the 
smooth  trajectory  rather  than  taking  the  abrupt  change  in  direction 
which  would  also  be  consistent  with  the  angular  momentum  of  the 
system.  That  is,  an  approaching  particle  cannot  suddenly  jump  to  a 
receding  trajectory.  The  finite  transition  may  be  viewed  as  a series 
of  infinitesimal  transitions  resulting  in  a gradual  distortion  of  the 
trajectory,  during  wliich  the  distinction  between  approaching  and  re- 
ceding paths  remains  clear  at  every  step.  For  the  finite  transition 

ij 

this  simply  implies  that  if  T at  is  an  approaching  path,  then 

Icl 

T is  to  be  chosen  as  the  approaching  path,  and  similarly  for 
receding  paths. 


I 


3.  Orientation  of  the  plane  F: 


By  th3  definition  of  the  angular  momentum,  the  plane  F in 
which  the  new  trajectory  lies  must  now  be  so  oriented  £S  to  be  per  - 
pendicular  to  the  angular  momentum  vector  of  the 

new  trajectory.  A difficulty  arises,  however,  due  to  the  fact  that 
it  is  not  possible  to  construct  a trajectory  through  B which  has  an 
arbitrary  angular  momentum  with  respect  to  the  point  0.  That  is, 


the  expression 


cannot  be  solved  for  except  in  the  rather  special  case 

that  is  perpendicular  to  . This  constitutes  a strong 

AO 

selection  rule  on  classical  transitions  in  which  the  internal  angular 

momentum  goes  from  ^ and  ^ to  and  ^ , 

the 

Unless  the  net  change  of^component  of  the  internal  angular  momenta 
along  the  radius  vector  is  zero,  the  transition  is  forbidden. 

When  the  transition  is  defined  quantum  mechanically,  however, 
the  angular  momenta  of  the  initial  and  final  states  are  not  completely 
specified,  so  that  in  the  correspondence  limit  a single  quantum 
niechar.’ical  transition  corresponds  to  a whole  family  of  classical 
transitions,  and  the  significance  of  this  selection  rule  is  somewhat 
modified.  This  may  be  seen  by  considering  the  two  general  cases 

r,Kf  ^ 

which  crise  when  we  solve  Eq.  2.2^12  for  CRab,)  . 

The  simplest  solution  of  this  equation  arises  when  we  consider 
a transition  in  which  the  angular  momentum  associated  with  the 
internal  motions  remains  unchanged,  as  is  the  case  for  example 
in  a pure  vibrational  transition.  In  this  case  Eq.  2.  2-12  immediately 
gives  the  result  that  is  coline ar  with  so  that 

the  plane  F coincides  with  the  xy  plane  and  the  angle  is 

zero.  Such  a transition  is  always  allowed. 


T 
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When  the  internal  transition  does  involve  a change  in  the  internal 

ang\ilar  momenta,  as  in  electronic  or  rotational  transitions,  it  is 

necessary  to  consider  the  internal  angular  momentum  of  the  quantum 

12 

mechanical  systems  in  the  light  of  the  correspondence  principle 

Since  the  quantum  mechanical  operators  toir  the  x,  y,  and  z 

13 

components  of  the  angular  momentum  do  not  commute  , it  is  not 
possible  in  quantum  mechanics  to  specify  the  angular  momentum  of 
a system  precisely . It  is  possible,  however,  to  specify  the  square 
of  the  total  angular  momentum  and  an  arbitrary  component,  say  in 
the  z direction. 

The  situation  is  illustrated  in  Fig.  2.  2-4. 


Fig.  2.  2-4  Classical  and  quantum  mechanical 
descriptions  of  the  angular  momentum. 


12.  L>.  1.  Schiff,  loc.  cit.  , p.  16. 

13.  H.  Eyring,  J.  Walter  and  G.  Kimball,  Quantum  Chemistry, 
Wiley,  1944.  p.  39. 
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Classically  the  vector  ,£^  is  precisely  specified  by  I , 

the  polar  angle  0^  and  the  azimuthal  angle  ■ Quantum 

mechanically  we  may  specify  only  the  length  I and  one 

component,  say  the  z compf^nent  given  by  the  angle  QL  ■ In 
t*he  correspondence  limit  the  azimuthal  angle  is  completely 

uncertain,  so  there  is  equal  probability  that  the  angular  momentum 
vector  lies  anywhere  in  the  cone  of  half  angle  Bo'  Thus  the 

single  quantum  mechanical  state  corresponds  to  a whole  family  of 
classical  states,  and  a quantum  mechanical  transition  of  the  angular 
momentum  vector  between  iwo  such  cones  corresponds  to  a doubly 
infinite  set  of  classical  transitions.  It  is  the  composite  behavior 
of  this  family  of  classical  transitions  which  corresponds  to  behavior 
of  the  single  quantum  mechanical  transition. 

Thus  if  the  quantum  mechanical  transition  in  which 

involves  a change  in  angular  momentum,  Eq.  2.2-12 

may  not  be  solved  uniquely  for  (2«b)  . Instead  there  is 

obtained  a family  of  solutions  corresponding  to  an  equal  probability 

distribution  of  the  classical  variables  cp^ , ’ft  , and  . 

All  members  of  this  family  which  are  perpendicular  to  R*  . . hence 

ao 

which  are  consistent  with  the  conservation  of  angular  momentum,  are 

kl 

used  to  construct  final  trajectories  T . Originating  from  every 

ij  kl 

point  of  T , then,  we  expect  a family  of  final  trajectories  T 

Illustration  of  this  principle  in  the  general  ca.se  is  difficult, 
because  the  coordinate  system  in  which  we  analyze  the  transition 
(Fig.  2.  2-1)  is  not  in  general  fixed  in  space  so  that  analysis  of  the 
angular  momentum  vectors  in  it  would  be  very  complicated.  Consider 
the  special  case  which  arises  when  ^ is  spherically  symme- 

tric. Then  the  entire  trajectory  T^"^  lies  in  the  xy  plane  and  the 
coordinate  system  is  fixed  in  space  so  that  the  z direction  may  be 
taken  as  the  axis  along  which  components  of  all  internal  angular 
momenta  are  known  precisely.  In  Chapter  III  it  is  shown  that 
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the  relative  motioa  of  rigid,  linear  polar  molecules  may  be  approxi- 
mated by  such  a spherical  potential,  so  the  case  is  of  considerable 
physical  interest. 

The  restriction  that  perpendicular  to  may 

be  expressed  in  the  form 


(2.  2-lbl 


where,  for  example 

This  result  may  be  derived  immediately  from  the  conservation  of  the 
component  of  the  angular  momentum  along  OB,  noting  that  the  trajec- 
tory angular  momenta  have  zero  component  along  this  direction. 
Equation  2.  2-16  implies  a relation  among  the  four  azimuthal  angles 
j 3.nd  , ®o  may  be  looked  upon  as  a third 

selection  rule  which  modifies  application  of  the  correspondence 
principle . 

kl 

That  this  selection  rule  limits  the  family  of  T diverging 
from  point  b,  but  never  entirely  forbids  the  transition  may  be  seen 
by  rewriting  Eq.  2.2-12  in  the  form 

1 iaT  1 t-  ] 


r ^ 

4-  I ib  I L \ 

— I JJi\  /3xWv  -t-  ~t<X^  \ 

- u:  I L ^ ^ ^ 


Clearly 

cf>^  = ^ cp^  = Co=t  fb  ~ \ 


Jk 


I 


44. 


is  alwa/s  an  acceptable  solution  .-if  this  equation,  for  arbitrary  , 

Thus  for  any  internal  transition,  one  or  more  solutions  of  the 
angular  momentum  equation  (Eq.  2.2*12)  may  be  obtained  which 
give  values  of  ./f satisfactory  for  determining  p CRJl) 
and  the  angle  ^ . If  there  is  only  one  resulting  trajectory,  its 
probability  is  given  by  the  change  in  the  transition  probability 
over  that  increment  of  the  initial  trajectory 
at  which  the  transition  occurred,  as  discussed  in  Sec.  2.  1.  If 
there  is  a family  of  resulting  trajectories,  the  probability  of  the 
entire  family  is  determined  by  this  change  of  the  transition  probabil- 
ity, and  the  distribution  of  probability  within  the  family  is  obtained 
by  assuming  that  all  values  of  the  azimuthal  angles  ^a.  / 
and  which  are  consistent  with  Eq.  2.2-12  are  equally 

probable. 

Thus  we  have  a complete  and  unique  method  for  constructing  all 
possible  trajectories  in  molecular  collisions  and  of  assigning  a 
relative  probability  to  each  of  them. 

In  the  general  case  of  transitions  involving  a change  in  the  internal 
angular  momentum  there  will  always  be  some  relation  of  the  oort 
given  in  Eq.  2.2-16  relating  the  geometric  location  of  the  transition 
on  the  trajectory  and  the  four  azimuthal  angles  , 9^ 

and  . Thus  if  we  ask  about  a particular  quantum  mechan- 

ical transition  at  some  specified  small  interval  along  the  trajectory, 
any  thuce  of  these  azimuthal  angles  may  be  considered  as  independent 
variables  say  and  9'“^  * For  each  choice  of 

these  three  angles,  we  may  solve  Eq.  2.  2-15  for  a value  of 

^ (.  , "-f'i  / ) and  determine  the  angle 

by  requiring  that  the  plane  F be  perpendicular  to  the  new  trajectory 
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angular  momentxim.  The  family  of  final  trajectories  resulting  from 
a transition  which  takes  place  between  time  t and  t + dt  along  the 
trajectory  may  be  thought  of  as  a function  of  the  three  azimuthal 
angles  and  the  time,  so  we  may  write 


According  to  the  correspondence  principle  there  is  equal  probability 
that  these  angles  have  any  value  between  0 and  2 77  , so  if  we 
wish  to  consider  any  function  of  the  trajectory^  ^ its  value 

averaged  over  the  family  of  trajectories  resulting  from  the  quantum 


transition  at  time  t may  be  written 


pLrr*)}  - ^77’ 


!■’> 


J 

o o « 


(2.  2-18) 


(2.2-19) 


For  example,  the  average  value  of  f'c*)  is  given  by 

^ ‘ ill 

o „ o 

The  probability  that  the  transition  occurred  in  the  time  interval 
d (i) 

dt  ii  given  by  ■ so  the  average  value  of  the  quantity 

Fcry  for  all  possible  final  trajectories  resulting  from  a single 
encounter  is  given  by 


F 


£p 

rrr 


P<h  1 I I 

-aaoo'o 


R-r^' 


<}l-& 


From  the  selection  rules  resulting  from  the  conservation  of  energy 

and  the  magnitude  of  the  anpi-'ar  momentum,  it  is  to  be  remembered 

kl 

that  those  final  trajectories  T are  forbidden  which  give  imaginary 
velocities  or  values  o£  ^ ~ whic-h  are  less  than  the 

inter  molecular  separation  K ♦ 
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m.  SAMPLE  CALCULATIONS  FOR  RIGID  LINEAR 
POLAR  MOLECULES 


In  this  chapter  we  illustrate  the  general  theory  derived  in 
Chapter  II  by  making  sample  calculations  relevant  to  the  study  of 
T OtcL  blO  nal  energy  transfer  in  collisions  between  rigid  linear  mole- 
cules containing  ideal  electric  dipoles.  In  order  to  be  specific, 
actual  numerical  results  will  be  given,  using  molecular  constants 
for  HCl  molecules. 

An  approximate  spherically  symmetric  potential  energy  function 
for  describing  the  relative  motion  of  the  molecules  is  presented  in 
terms  of  London's  solution  of  the  adiabatic  clamped  equations.  The 
form  of  these  potentials  suggests  why  the  viscosity  of  polar  gases 
actually  can  be  described  by  the  present  theory  for  spherical 
molecules . 

The  various  transition  probabilities 
and  " ^ ^ introduced  in  Chapter  II  are  calculated  and  discussed 

for  the  special  case  of  constant  velocity  linear  trajectories.  Inves- 

d.  f ~ 

tigation  of  ^ L ’ indicates  that  the  distance  of  closest 

approach  in  an  encounter  is  an  unlikely  place  for  the  inelastic  trans- 
fer to  occur. 


3.  1 The  Isolated  Molecule  Wave  Functions  and  the  Inter  molecular 
Potential. 

a)  The  wave  functions. 

Consider  a diatomic  molecule  with  atoms  of  mass  mj^  and  m, 
separated  by  a distance  d . Suppose  that  at  the  center  of  mass  of 
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the  molecule  there  is  an  ideal  electric  dipole  p oriented  parallel 
to  the  molecular  axis. 

If  d is  held  fixed,  the  solutions  of  Eq.  2.  1-2  describing  the 
motion  of  the  molecule  relative  to  its  own  mass  center  are  simply 
the  rigid  rotator  wave  functions  and  energy  levels \ Thus  for 


molecule  a, 


(3.  1-1) 


Ei  " 


(3.  1-2) 


and  similarly  for  molecule  b.  Here  _I  is  the  moment  of  inertia, 
which  in  terms  of  the  reduced  mass  /i  is  given  by 


X ^ 


/TV,  r/Wv 


(3.  1-3) 


w'”'*-  2 

and  the  i C^,<p*.)  are  the  spherical  harmonics  defined  as 

The  Kj  are  the  associated  Legendre  functions,  give  >'a  by 


j 1^1 


I *■  j^i 


( 1-7-^)^ 


b)  The  potential  enerev  of  interaction. 

The  energy  of  interaction  between  two  ideal  dipoles  of  strength 


1.  For  an  example  of  the  quantum  mechanical  solution  of  the 
motion  of  a rigid  rotator,  see  for  example,  H.  Eyring, 

J.  Walter,  and  G.  Kimball,  loc.  cit.  , p . 72. 


2.  This  definition  is  taken  from  MTGL,  p,  905. 
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and  is  given  in  vector  fornn  by 

(h.  ’fo  - ‘tj.  /. 


(3.  1-4) 


^ iRaki-  " 

In  a coordinate  system  in  which  the  line  connecting  the  mass  centers 
of  the  two  molecules  is  taken  as  the  polar  axis,  so  that  6>o,  and  &b 
are  polar  angles  giving  the  inclination  of  molecules  a and  b relative 
to  this  axis,  and  and  <^c  the  azimuthal  angles  about  this 

line,  takes  the  special  form 

^ — (3.  l-4b) 


RaJ- 


The  total  energy  of  interaction  between  polar  molecules  may  be 

4 

described  approximately  by  the  Stockmayer  potential 


=.  yeL(i£J'  ' '|£ 


(3.  1-5) 


i- 


>£t 


r< 


Ui?  I 


Kat 


I “ 


The  first  brrcketed  term  in  this  potential  is  a Lennard -Jones  type 
spherically  symmetric  potential  describing  the  hard  core  repulsive 
energy  and  the  attractive  dispersion  energy.  The  remainder  is  just 
the  dipole -dipole  interaction  energy.  The  constants  (T"  and  6- 
are  characteristic  of  each  molecule  and  are  determined  experimen- 
tally . 

The  molecular  constants  for  HCl  which  will  be  used  in  sub- 

5 

sequent  calculations  are  as  follows  : 


3.  MTGL,  p.  851. 

4.  W.  H.  Stockmayer,  J.  Chem.  Phys.  2.>  398,  (1941). 

5.  The  values  of  I and  M.  are  taken  from  G Kerzberg, 

Spectra  of  Diatomic  Molecules,  Van  Nostrand,  1950,  p.  534. 
Those  of  [%.]  from  MTGL,  p.  1112. 
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I = 2.6422 

X 

10-"° 

2 

gm.  cm. 

fi  = 1.62645 

/ 

10-"" 

gm. 

p = 1.034 

X 

o 

1 

00 

e . s . u. 

= 3.305 

X 

10 -« 

cm. 

(3.  1-6) 


= 360  ° |< 


Here  are  constants  obtained  by  fitting 

viscosity  data  to  the  Lennard-Jones  potential,  and  would  not  be 
expected  to  be  precisely  the  same  as  the  constants  obtained  by  fitting 
equation  of  state  data  to  the  Stockmayer  potential.  Since  Stockmayer 
constants  are  not  available  for  HCl  the  L.  J.  values  will  be  used 
as  a rough  estimate  of  the  spherically  symmetric  potential. 


3.  2 The  Effective  Potential 

It  is  not  possible  to  solve  Equ.  Z.  1-3  exactly  for  the  adiabatic 
clamped  wave  functions  4^^  and  the  effective  potential  7 ( 1^*’)  • 

London^,  however,  has  obtained  approximate  solutions  of  the  problem 
which  allow  us  to  construct  approximate  effective  potentials  which 
depend  only  on  the  distance  of  separation  | I . London  used 
a modified  perturbation  method  which  he  called  the  method  of  unsharp 
resonance,  but  the  solution  is  equally  well  viewed  as  a variation 
method. 

When  the  mole  cules  are  far  separated  the  wave  fvinction  of  the 


6.  F.  London,  Z.  Physik,  63,  Z45  (1930). 
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system  may  be  written  as  a product  of  spherical  harmonics 


= U^o  ^ (3.2-1) 


For  convenience  we  write  this  function  for  the  initial  state  a 
where  the  superscript  denotes  that  it  is  an  unperturbed  function, 
and  the  subscript  stands  for  the  four  quantum  numbers  j , m , 

and  . For  a trial  function,  London  forms  a linear  combina- 
tion of  ~]po  and  the  twelve  (at  most)  unperturbed  functions 

with  which  it  combines  in  the  sense  that 


( ii)°  I (P“)  o 


These  twelve  functions  are 

v!"  - Ob*-', 

Ipt  - J -'i 

4Jj  - U'’(]a-l,/Wa+' ; ib-', 

^ '^a'l  y }b'' / 


^7  --  ^b) 

144-  = •>) 

iPZ  = ^ b->-> , 'rr>  t) 

l/^/V  " r^b.-ry 


(3.  2-2) 


Functions  with  negative  j values  or  values  of  | m|  greater 
than  the  corresponding  j are  taken  to  be  zero. 

The  best  value  of  the  energy  obtainable  from  this  linear  combin- 
ation is  given  by  solution  of  the  13 'th  order  secular  equation 

I ( ) ■*■  C Ec  ~^o  " (3.2-3) 


! 


^ ^ O,  i, 


/ a. 


Here  the  Ej^  are  the  energies  corresponding  to  the  unperturbed 
wave  functions  . If  relative  coordinates  are  chosen  so  that 


the  intermolecuiar  axis  is  taken  as  the  z direction,  the  secular 
determinant  takes  on  a relatively  simple  form.  All  matrix  elements 


vanish  except  those  in  the  first  row  and  the  first  column,  and  the 
diagonal  elements  become  equal  in  groups  of  three.  London  shows 
that  the  secular  equation  may  be  reduced  to  the  fifth  order  equation 


I ^ a; 

f' 


f-t.  9-^ 


(3.2-4) 


Here  the  6:,;  are  given  by 


(3.2-5) 


L - 4^  ■? 


, = i.'-f  -h  + 

tv  ^ L ^ J."t, J 


^ evil  ^^1  + C I ^ 

‘ ^ (3. 2-6) 

The  are  given  explicitly  in  terms  of  the  quantum  numbers  in 

•7 

MTGL  *.  For  example,  the  coefficient  Ctf  is  given  by 

V [ ][(<}!.  '7 

■h±  I 1-nria.)  ( ^oi-nrr^~l)  ( ^ fc  r/V7i^  ,W^  ( /b 

[>'3^  -'IL  VQb  + i}‘J 


The  solutions  of  Eq.  3.  2-3  fall  into  two  general  classes,  categor- 
ized according  to  their  behavior  for  large  values  of  R , . For  non 

cLD 

resonance  interactions  ( j ^ j + 1)  the  physically  significant 

3.  D ^ 

solution  is  negative  and  is  proportional  to  1 • It  is  shown  in 

MTGL  that  the  coefficient  of  this  term  is  given  by 


7.  MTGL  p.  1000. 
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( Rah  ~>ao) 

Vyufn  rj^ 


(3.2-7) 


For  resonant  interactions  (j  = j.  + 1 , and  p = p,  = p) 

a D a b 

there  are  two  solutions  of  Eq.  3.2-3  which  for  large  differ 

only  in  sign  and  are  proportional  to  1/r\.  In  MTGL  the 

ab 

coefficients  are  given  as 


(3.  2-8) 


Solutions  of  Eq.  3.  2-3  for  these  two  cases  are  tabulated  in  Table 

3.  2-1  and  are  plotted  in  Fig.  3.  2-1,  for  the  interaction  of  two  HCl 

nnolecules.  Curve  A gives  the  non  resonant  interaction  when  both 

molecules  are  in  their  ground  states.  The  pair  of  curves  labeled  B 

give  the  interaction  energy  for  the  resonant  case  j =1,  m =0; 

a a 

= 2,  = 1.  The  pair  labeled  C are  for  the  resonant  case  j =2, 

= 0;  For  comparison,  curve  D gives  the 

attractive  part  of  the  Leonard -Jones  potential,  which  is  proportional 

6 

to  1/R^j^.  Molecular  constants  are  taken  from  Eq.  3.  1-6. 
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10.5  -0.0685  -0.1^3  ''10.0  4.432  3. 581 

11.0  -0.0520  -0.146  -11.0  4.306  3.474 

11.5  -0.0393  -0.112  -12.0  4.194  3.375 

12.0  -0.0310  -0.087 
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a)  Non  Resonant  Interactions 

The  non  resonant  interactions  are  best  examined  by  looking  at 

the  coefficient  of  the  I /R  term  which  describes  the  interaction 

ab 

at  large  separations.  The  constant  for  several  sets  of 

- 8 

quantum  numbers,  and  for  R measured  in  units  of  10  cm.  , 
is  given  in  Table  3.2-2. 


Table  3.  2-2 

Coefficients  of  the  interaction  energy  term  proportional 

to  1 /R  , in  non  resonant  collisions 
ab 


Internal  states 
= 0 , m = 0;  j r 0, 

a b 

= 5.  m^  = 2;  = 5, 

= 10,  m^  = 0;  - 10, 

= 5,  ma  = 0;  = 10, 


[erg  % 

o 

II 

9.06  ^ 10'^^ 

u 

-10 

m^  = 2 

1.  12  ^ 10 

-10 

™b 

1. 34  A 10 

o 

II 

6.04  / 10*^^ 

This  may  be  compared  to  the  coefficient  of  the  term  in  the  Lennard- 

Jones  potential  which  varies  as  1/R  , , 

ab 

B,  - ^ ^ 

For  the  values  given  in  Eq.  3.  1-6  for  HCl,  ~ 

Since  B represents  some  sort  of  an  average  of  the  dispersion 
energy  and  the  various  dipole -dipole  forces^  it  appears  iiKeiy  mat 
the  dispersion  energy  is  at  least  oi  the  same  order  of  magnitude  as 
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the  various  non  resonance  dipole -dipole  forces.  London  investi- 
gates this  question  in  detail  for  several  polar  gases,  and  concludes 
that  none  of  the  forces  between  HCl  molecules  which  arise  directly 
from  the  electric  dipoles  are  as  large  as  the  dispersion  force. 

b)  Resonant  interactions. 

The  resonant  interactions  are  of  particular  interest,  and  are 
best  investigated  by  looking  at  the  coefficient  a^  (Eq.  3.  2-6b). 

First,  however,  note  by  comparison  with  Eq.  3.  l-4h  for  the  static 
interaction  of  two  dipoles,  that  a^  represents  some  sort  of  an 
average  over  the  angular  factor 

iSs  & >-  &t, 

which  gives  the  orientation  of  the  dipoles  relative  to  the  intermolecular 
axis.  The  absolute  value  of  this  factor  varies  from  a maximum  of  two 
when  the  molecules  are  parallel  ( 6>a  - 0,71  \ 9b  - O , to  a minimum 

of  zero  when  they  are  perpendicular  ( ■-  - ^)<  Hence  we  ex- 

pect a^  to  lie  somewhere  in  this  range  . For  some  typical  resonant 
interactions  a^  has  the  values  given  in  Table  3.  2-3. 


I 

I 

i 

I 8.  F.  London,  Trans.  Faraday  Soc . , 189  , 8 (1937) 


I 


'iiiiii' . fi. 
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Table  3.2-3 

London  Factor  for  the  Coefficient  of  the  Interaction 
Energy  Term  Proportional  to  l/R^^^  in  resonant 
Collisions. 

Internal  states  a^ 


^a 

= 1. 

m 

a 

= 0;  j^ 

= 2, 

^-b 

^ 1 

— A 

0.  586 

■^a 

= 1. 

m 

a 

==  0;  jb 

= 2, 

= 2 

0.436 

ja 

= 1, 

m 

a 

O 

II 

= 2, 

= -1 

0.  169 

ja 

= 4, 

m 

a 

II 

o 

u... 

cr 

= 5, 

= 0 

0.  531 

r Apparently  the  averaged  dipole-dipole  energy  is  of  the  order  of  25% 

or  less  of  the  maximum  interaction  energy  of  the  static  dipoles. 

9 

Krieger  performed  a very  interesting  calculation  of  the  viscosity 
of  polar  gases,  by  using  a sphericaEy  symmetric  potential  function 
related  to  the  Stockmayer  potential.  He  assumed  that  the  energy  of 
interaction  between  polar  molecules  could  be  expressed  in  the  form 


ye 


/12T  i 


Here  p is  the  experimental  value  of  the  dipole  moment,  and  ^ 
and  e are  adjustable  constants  used  to  fit  the  experimental  data. 
This  potential  energy  function  corresponds  to  the  Stockmayer  poten- 
tial when  the  dipoles  are  oriented  in  the  parallel  orientation  of  minimum 
energy,  as  in  the  accompanying  diagram. 


9.  F.  J.  Krieger,  Project  Rand  Repox-t  RM-646  (1951) . 
See  also  MTGL  p.  597. 
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R., 

For  this  orient**tion  the  angular  factor  of  the  Stockmayer  potential 
(Eq.  3.  l~4b)  is  equal  to  minus  two. 

Although  this  theory  ignores  the  effects  of  the  internal  degrees 
of  freedom,  Krieger  is  able  to  fit  the  experimental  viscosity  data 
of  polar  gases  over  the  entire  experimental  temperature  range  to 
within  experimental  error.  This  surprising  success  might  be  inter- 
preted as  follows: 

1.  An  appreciable  fraction  of  the  collisions  in  a gas  of  rigid 
rotators  are  resonant  collisions.  ( 3*,  - - >)  The  fraction 

IS  given  by 


r 


of  resonant  collisions,  Tres 

r~ 





J 


(3.2-9) 


e 


[ 2 e 

L 3 


J 


-»rA.r 

For  large  values  of  I and  T this  may  be  evaluated  by  replacing 
the  sums  by  integrals.  In  this  case 


■f  r«s 


. 1; 


jr 

'/  s. 


(3.  2~9a) 


For  collisions  between  HCl  molecules  at  300  K.,  this  gives 

4-  t-ei  - - -2-  S 3 

2. In  resonant  collisions  the  translational  energy  is  conserved, 
so  that  the  collision  is  sort  of  a quasi  clastic  collision.  In  Eq.  1. 2-2 
for  the  Uhlenbeck  iorm  of  the  viscosity  coefficient  the  term  multiplied 
by  the  net  change  of  internal  energy 


vanishes,  and  the 
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property  we  average  over  all  possible  encounters  is  just 


— CO-4-' 


yxj 


as  in  the  expression  far  spherical  molecules,  Eq.  1.2-1. 


3.  Apparently  the  quasi  elastic  resonant  encounters  are  of  such 

importance  in  viscosity  that  Krieger  was  able  to  predict  the  esBcatial 

behavior  of  the  gas  simply  by  including  the  long  range  energy  term 

3 

proportional  to  have  seen  that  in  all  non  resonant  inter- 

actions the  coefficient  of  this  term  should  be  zero,  and  that  in 
resonant  collisions  it  should  depend  an  the  quantum  states  of  the 
colliding  molecules.  In  resonant  collisions,  moreover,  it  shoiild 
be  considerably  smaller  than  Krieger 's  value  of  2"P  , and  should 
be  repulsive  in  some  collisions,  and  attractive  in  others.  Adjustment 
of  the  coefficients  ^ and  €:  in  Krieger 's  theory  is  able  to 

compensate  for  all  of  these  approximations. 

Although  these  considerations  have  been  in  terms  of  our  study 
of  linear  molecules,  whereas  Krieger’s  measurements  included 

non  linear  molecules,  the  general  behavior  would  be  expected  to 

10 

be  the  same.  Margenau  and  Warren  have  investigated  the  first 

order  interactions  of  ideal  dipoles  imbedded  in  symmetric  tops,  and 

3 

find  interactions  proportional  even  in  the  non  resonant  cases. 


c)  Suggested  procedure  for  considering  the  viscosity  of  "olar  gases 
From  the  considerations  in  this  section,  of  the  effective  inter- 
molecular  potential  arising  between  linear  molecules  containing 
electric  dipoles,  we  would  suggest  the  following  procedure  foi  pre- 
dicting the  viscosity  of  gases  composed  of  such  molecules: 


10.  H.  Margenau  and  D.  T.  Warren,  Phys.  Rev,,  51,  748  (1937) 
See  also  ^4TGI->,.  p.  1004. 
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Consider  the  gas  as  a mixture,  with  each  set  rf  rotational  quantum 
numbers  defining  a separate  molecular  species.  The  fraction  of  the 
molecules  which  are  of  species  specified  by  the  rotational 

quantum  numbers  ja,  and  mo,  is  given  at  thermal  equil- 
ibrium by  ^ 


- SJUT 


(3. 2-10) 


When  two  molecules  of  species  a and  b collide  such  that 

the  dynamics  of  the  encounter  are  determined  by  the  inter- 


molecular  potential  function. 


cp  = Y 


- (£)-] 


(3.2-11) 


Where  <T~  and  6 are  experimental  constants  determined 
by  fitting  experimental  values  of  the  second  virial  coefficient  as  a 
function  of  temperature  to  the  Stockmayer  potential  (Eq.  3.  1-5). 
The  fraction  ot  the  collisions  in  the  gas  which  are  governed  by  this 
potential  is  (^1  where  ■fre.>  is  given  by  Eq.  3.  2-9a  . 

When  species  a and  b are  such  that  j = j,  + 1,  the 
potential  energy  of  interaction  is  taken  to  be  of  form 


'p-  "4(€r-(£r] " 


(3.  2-12) 


here  p is  the  experimental  value  of  ihc  dipole  moment,  and  the 
quantity  ^ is  related  to  the  quantities  a^  defined  in 
Eq.  3.2-6b. 

L>et  us  call  a collision  in  which  either  of  the  molecules  is  in 


3ta.te  a and  the  other  in  state  b an  a-b  collision. 


If  we 


label  our  molecules  such  that  j = j + 1,  then  in  half  of  the 

SL  D 

a-b  collisions  we  take 


1 
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i 


S 

f 


and  in  the  other  half 

f(  jc, 

The  fraction  of  the  collisions  in  the  gas  which  are  a-b  collisions 

is  just  2f  f , where  the  factor  of  two  results  from  the  fact  that 
3i  D 

either  molecule  may  be  in  the  higher  energy  state. 

The  dynamics  of  encounters  governed  by  the  .Lennard- 
Jones  type  potential  of  Eq.  3.2-11  have  been  calculated  and 
tabulated.  The  dynamics  of  collisions  governed  by  the  potential 

g 

of  Eq.  3.2-12  are  given  in  Krieger's  v;ork  , for  negative  values 

of  f , If  K rieger's  calculations  were  extended  to  include  positive 

values  of  ^ it  would  not  be  a difficult  matter  to  test  this 

theory  of  the  viscosity  of  linear  polar  molecules. 

The  viscosity  of  synametric  top  molecules  with  imbedded  electric 

dipoles  could  be  calculated  in  a manner  very  similar  to  this,  using 

9 

the  effective  potentials  calculated  by  Margenau  and  Warren  . The 
chief  difference  wculd  be  that  a larger  fraction  of  the  encounters 
would  be  governed  by  potentials  which  contained  a term  proportional 

3. 3 The  Probability  of  Transition 

, , K . K , ^ . 

The  three  transition  probabilities  introduced  in  Chapter  II  were 


11.  MTGL  p.  1132. 
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(2.2-6) 


Here  P°6S  was  derived  from  a time  dependent  perturbation 

procedure  in  which  the  unperturbed  functions  were  taken  to  be  simple 
products  of  the  Isolated  molecule  wave  functions  JJ.L('i)  and  M-b  <■0-)  . 

vrAB  obtained  by  using  the  adiabatic  clamped  functions 
as  a time  dependent  unperturbed  set.  The  function 
suggested  in  Sec.  2.2  as  a senoi-classical  analogue 

of  K^'Ut) 


In  order  to  better  understand  these  three  functions,  we  calculate 

each  of  them  for  the  special  case  in  which  the  equivalent  one  particle 

trajectory  is  a straight  line  traversed  at  constant  velocity  (Fig.  3.  3-1) 

Although  this  approximation  would  not  in  general  be  valid  if  we  were 

seeking  precise  numerical  results,  it  will  be  adequate  for  illus»^rating 

and  comparing  the  general  behavior  of  these  functions.  The  advantage 

of  the  approximation  is  that  the  time  dependence  of  the  matrix 

elements  of  ^,u-  j originating  in  the  factor  1 /R  ^ now 

~V\- 

appears  simply  as  L -i 

For  the  purpose  of  these  approximate  calcuJ.ations  we  also 
replace  the  adiabatic  clamped  functions 


and 
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Fig.  3.3-1  The  simplest  case  of  an 
equivalent  one  particle  trajectory. 


— <•  S . • 

energies  t C^at,)  by  their  zero  order  approximations  Mb 

and  E ^ . 

With  these  approximations,  the  transition  probabilities  may  be 


written 


(■t) 

t 4 


r — ^ 


L . ] 

“i; 


p'y  = [3C  f ■^a’/ 


T.  T ^ 


64. 


R . ^ 

1 it 


HI  ct  2 ^0/ 


(3.3-3) 


The  various  q\iantities  appearing  here  are  as  fellows; 
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^ f/  - - e. 


2 - 
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^ ^ j 


'-  ^ _ . i‘  e / po. 

' ' *■  ~ ^ ^E"'' 


<■1 


e ( Po.  ipt,  )(  3 
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Ka 


^ M-  to  j ( Mo-  M b ^ 


fo.  ft, 


~ i2  Cog,  c-v7^  & 


•4ijj  W'^0 

>J  r,,(ft.%;  Ij/o.,f0 


o o ^ 

Here  the  quantum  numbers  abbreviated  by  ij,  and  kl  are 

i ^ ^(Ky  t<  ^o., 

3 3 t) , /W  b £ "^  >•  b , /■'-^ 


We  also  use  the  result  that 

<x>  - o H 

r ^ 


12 


r i y-t,  <7 


- 77  I , , 

z:  W ^ c.; 


" ' 


— vA/' 

|(0 


is  the  first  order  Hankel  function  of  the  first 
13  . 


Where  H; 

kind,  tabulated  in  Jahnke  and  Emde^’’.  The  integrals  <£^c> 


12.  G.N.  Watson,  A Treatise  on  the  Theory  of  Bessel  Functions 
Camb.  Univ.  Press,  Seconded,  p.  172 

13.  E.  Jahrke  and  F.  Emde,  fables  of  Functions,  Dover,  1945,  p.  133. 
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may  be  calculated  directly  from  the  properties  of  spherical  harmonics 

14 

or  they  may  be  found  listed  in  convenient  form  by  Margenau 

In  Table  3.  3-1  and  Fig.  3.3-2  we  give  calculated  values  of 
the  three  transition  probabilities  for  the  case  when  the  initial  relative 
velocity  g (taken  tc  be  constant  throughout  the  encounter),  the 
collision  parameter  b,  and  the  change  in  internal  quantum  sf.ates 
are 

5 

g = 3.  1924  X 10  cm.  /sec. 


b 

=4.89  X 

10 

-8 

cm. 

ja 

= 4,  m = 
a 

0; 

ib  = 

= 0 

A a 

= 5,  /<a  = 

0; 

Xb  = 

% 

= 0 

The  constants  correspond  to  a typical  encounter  for  HCl  molecules 
at  room  temperature. 

d''' 

The  principal  deviations  of  C from  the 

general  form  of  r are  not  too  difficult  to  understand.  The 

i-l  -^Kt 

oscillatory  motion  of  these  two  probabilities,  with  period  equal  to 
i'./ EiV  is  21  typical  quantum  mechanical  phenomenon  which 

may  be  interpreted  either  as  a diffraction  effect  resulting  from  the 
wave -like  nature  of  matter,  or  as  a consequence  of  the  uncertainty 
principle.  From  the  point  of  view  of  the  uncertainty  principle,  we 
are  not  allowed  to  specify  the  time  interval  during  which 

kJ 

the  energy  cliar.j^e  took  place  more  precisely  than  about 


14.  K.  Margenau,  F.ev.  Mod.  Phys.  11, 


1 (1939). 
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Table  3-3-1 
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Traneition  i*robabllltlei  alon«  a Straight  Lins  Trajectory’ 
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0.0218 

0.0000 
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0.0072 
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0.0085 
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0.0000 

0.0000 

0.0116 

-5.0 

0.022b 

C.0005 

0.0195 

-4.5 

0.0296 
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0.0340 

-4.0 

0.0709 

0.0023 

0.0621 
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0.142 

0.0053 

0.0006 

0.1190 

-3.0 

0.205 

0.0127 

0.0024 

0.2452 
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0.497 

0.0322 
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0.5325 
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0.0950 

0.0084 

i>.221 
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0.0230 

2.867 
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5.97 

0.4214 
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-0.5 
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0.50iW 
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11.82 

0.0 

19.1 

0.3344 

0.1496 

14.95 

0.5 

16.6 

0.8667 

0.1860 

11.82 

1.0 

7.41 

1.276 

0.2423 

6.4o9 

1-5 

1.21 

0.6576 

0.2782 

2.867 

2.0 

1.16 

0.0950 

0.2877 

1.221 

2.5 

1.51 

0.268 

0.2961 

0.5325 

3.0 

0.6i4 

0.4292 

0.2967 

0.2452 

3*5 

0.0361 

0.3278 

0.2985 

0.1190 

4.0 

0.263 

0.2508 

0.0621 

^■5 

0.490 

0.2909 

0.0340 

5.0 

0.320 

0.3242 

0.0195 

5.5 

0.178 

0.3037 

0.0116 

6.0 

0.269 

0.2863 

0.0080 

6.5 

0.351 

0.2972 

0.0062 

7.0 

0.301 

0.306S 

0.0031 

7-5 

0.232 

0.3003 

0.0021 

8.0 

0.274 

0.2946 

0.0014 

00 

0.2991' 

0.2991 

0.2991 
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This  means  that  wc  are  not  allowed  to  ask  about  the  probability  of 
transition  during  time  intervcils  less  than  that,  corresponding  to  about 
one  cycle  of  the  curves  of  Fig.  3.  3 -1 

The  large  maximum  in  about  ■6-0  may  be  understood 

by  noting  that  actually  gives  the  transition  probability  for 

an  experiment  in  which  the  c.  ’ision  is  allowed  to  proceed  until 
time  t,  that  the  perturbation  is  suddenly  turned  off.  The  sudden 
turning  off  of  the  perturbation  itself  induces  transitions,  corresponding 
to  the  fact  that  at  that  moment  the  true  wave  function  of  the  system 
has  large  contributions  of  other  unperturbed  states,  simply  because 
of  the  adiabatic  distortion.  For  comparison,  curve  D of  Fig.  3.3-2 
is  a plot  of 

u:  i I ^ 

I 

the  first  order  perturbation  estimate  of  the  contribution  of  the  kl 
state  to  the  wave  function  in  an  adiabatic  encounter,  This  curve 
accounts  very  well  for  the  wide  maximum  in  • 

The  transit.ion  probability  corresponds  more  nearly  to 

an  experiment  in  which  at  time  t the  relative  motion  of  the  molecules 
is  stopped  and  they  are  separated  adiabatically  to  infinity.  The  fact 
that  there  may  still  be  something  of  a maximum  in  the  vicinity  of 
the  point  of  clo.Test  approach  (t  = 0)  may  be  due  to  the  approxima- 
tion of  evaluating  the  energy  and  matrix  elements  in  terms  of  the 
The  fact  that  and  approach  the  same  value  at 

t = CO  is  understandable,  since  in  this  limit  the  adiabatic  functions 
become  identical  with  the  simple  product  functions.  Mathematically 
the  equivalence  rr.zy  be  demonstrated  by  integrating  Eq.  3. 3-1  by 
parts;  consistent  with  the  assumptions  we  have  made  in  the  evaluation 
of  P.  , it  fcllo'ws  immediatelv  that  , far  all 
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b)  Dependence  cf  the  transition  probability  on  g,  b,  and  ttj 

To  show  how  strongly  the  transition  probability  depends  on  the 
relative  velocity,  g,  the  collision  parametei;  b,  and  the  energy 

-'t  pt'i 

exchange  At",;)  , we  have  calculated  for  collisions  between 

HCl  molecules,  varying  each  of  these  quantities  separately.  The 
results  are  given  in  Table  3.  3-2. 

From  these  calculations  we  see  that  changing  each  of  these  quan- 
tities by  less  than  a factor  of  two  changes  the  total  probability  of 
transition  by  a factor  of  about  one  hundred.  In  an  actual  calculation 
of  the  cross  sections,  this  would  mean  that  the  main  contribution  to 
each  differential  cross  section  would  be  expected  to  come  from  a 
relatively  narrow  range  of  relative  velocities  and  collision  parameters. 

c)  Location  of  the  transition  along  the  trajectory. 

In  the  semi-classical  theory  of  coilisiuns  it  is  necessary  to  know 
where  along  the  trajectory  a given  transition  has  taken  place.  In  our 
semi-classical  description  of  collisions  this  is  given  by  the  shape  of 
T)  or  by  the  related  function  tf<r  L L . 

r~N  ^ 

I t:;  gives  the  probability  that  the  transition  has  occurred  up 

I - 

4 - 

to  time  t,  and  ' r >■'  ■ gives  the  probability  that  the  transition 

occurs  in  the  time  interval  about  t.  From  Eq.  3 3-3  it  may  be 

seen  that  to  the  approximations  used  in  deriving  that  equation,  each 
of  these  is  given  by  "shape  factors"  which  depend  only  on  the  reduced 
variable  , z = gt/b.  Thus  they  may  be  written 


/ it) 

i L j 


(3.3-4) 


and 
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where 


and 


z r 1^1  di 
^ ^ la,Ui- 


7.  JtL 


To  illustrate  the  behavior  of  these  functions,  we  give  calculated 
values  of  S(z)  in  Table  3.  3-3  and  plot  the  results  in  Fig.  3. 3-3.  The 
most  interesting  feature  of  this  curve  is  the  plateau  in  the  vicinity  of 
t = 0,  the  distance  of  closest  approach.  The  significance  of  this  plateau 
is  seen  even  more  clearly  in  Fig.  3.  3-4  where  we  plot  the  function 
S'(z).  From  this  curve  we  see  that  the  probability  that  the  transition 
occurs  during  an  infinitesimal  time  increment  about  t = 0 is  zero. 

The  maximum  probability  of  transition  occurs  when  . 

The  location  of  the  maximum  of  the  function  exactly 

b,  ;v1  is  of  course  due  to  our  approximate  form  of  P 
but  the  zero  at  the  distarc<*  of  closest  approach  appears  even  in  the 
accurate  quantum  mechanical  formula.  Thus 


> o «%»  *■  ■' 


/i,  .*<■  (4# ; 
V '+'.^  I dt  I 


Since  the  distance  of  closest  approach  is  a stationary  point  of  Rab* 

( I 3^1  — O at  t = 0.  Thus  the  probability  that  the 

transition  occurs  in  the  time  interval  dt  about  the  point  t=0  vanishes. 

One  must  be  careful  in  interpreting  this  result,  however,  not  to 
conclude  that  the  segment  of  the  trajecto~y  about  the  distance  of  clusest 
approach  is  always  an  unlikely  place  for  a transition  to  occur.  In 
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givet  the  prouaole  distribution  of  tr.nsitioas  during  the  time  of  the  collision. 


Sable  3.3-3 


Shape  Jactori  Aesoclatea  with  for  a Straight  Line 

Trajectory 


If 

(II 

S (s) 

S'(s) 

0.00 

0.500 

0.000 

0.05 

0.504 

u.  j.(3 

0.10 

0.517 

0.335 

0.475 

0.15 

0-537 

O.PQ 

0.5S4 

0.527 

0.25 

0.596 

0 . bo6 

0.^0 

0.630 

0.712 

o.lk) 

0.703 

O.7I8 

0.50 

0.771 

0.641 

0.60 

0.830 

0,526 

0.70 

0.876 

0.407 

O.BO 

0.511 

0.302 

0 on 

0.337 

0.218 

1.00 

0.956 

0.155 

1.10 

0.969 

0,109 

1.20 

0.978 

0.076 

0.984 

0.053 

1.40 

0.989 

0,037 

1.50 

0.992 

0.026 

1.60 

0.994 

0.018 

1.70 

0.996 

0.013 

1.80 

0.997 

0.009 

1.90 

0.998 

n AAT 

0 w 1 

2.00 

0.998 

C.OO5 

* Values  of  8 (*)  and  S'(*)  for  negatire  1 aay  be  obtained 
by  the  relations  S (-z)  = 1 - S (z)  and  SH-t)  = S (+z). 
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some  cases,  as  in  head  on  collisions,  the  particles  may  spend  a long 
time  in  the  configuration  of  closest  approach  and  the  probability  of 
transition  in  this  region  may  not  be  small.  This  may  be  seen  by  writing 
as  a function  of  the  arc  length  s.  Then 

^ . f ^ H/ 


AIC<?.1 


d 

dS 


lF%>! 


(3.  3-7) 


Now  in  trajectories  which  are  nearly  like  the  straight  line  constant 

■ ^ 

velocity  paths  used  in  this  section,  is  indeed  zero  at  the  distance 

diS^ 

of  closest  approach,  so  that  the  probability  of  a transition  occurring 
in  the  element  of  trajectory  ds  about  the  distance  of  closest  approach 
is  indeed  zero,  as  pictured  in  Fig.  3.  3-4.  In  a head-on  collision, 
however,  the  arc  length  s becomes  identical  with  the  intermolecular 
distance  Rab«  ^.nd  the  function  is  now  large  in  the  vicinity  of  the 

distance  of  closest  approach.  (Actually  it  is  t.ol  defined  just  at  the 

distance  of  closest  approach  because  of  a cusp  in  the  function 

' d ^ 

at  that  point,  but  on  both  sides  ot  the  cusp  will  in  general  be 

large.)  Thus  in  the  limit  ot  a head  on  collision,  the  probability  of  a 
transition  in  the  element  ds  about  the  distance  of  closest  approach  may 
become  large.  It  is  interesting  to  note  that  in  the  limit  of  b = 0,  a 
head-on  collision,  the  two  maxima  of  Fig.  3.  3-4  move  toward  the 
origin,  indicating  the  likelihood  of  a transition  in  this  region. 

The  trajectory  along  which  we  are  carrying  out  the  time  dependent 
perturbation  calculation  is  symmetric,  about  t - 0,  so  in  some  sense 
it  is  true  that  the  average  location  of  the  transition  in  grazing  collisions 
is  at  the  distance  of  closest  approach.  In  kinetic  theory  applications, 
however,  we  are  interested  in  the  angle  of  deflection  resulting  from 
the  transition,  so  that  it  is  by  no  means  valid  to  expect  that  the 
distance  of  closest  approach  would  be  an  appropriate  average  location 
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for  all  transitions. 

d}  The  transition  probability  for  rtoonant  collisions 

The  considerations  of  the  probability  u£  trancition  thus  far  in  this 
section  have  been  limited  to  non  resonance  interactions,  since  we  have 
assumed  that  ^ io  all  of  our  formulae  and  calculations. 

15 

The  classical  analogue  of  a quantum  mechanical  system  in  resonance 
pictures  the  extra  quantum  of  energy  as  being  transferred  rapidly 
back  and  forth  between  the  colliding  molecules,  so  that  in  a measure- 
ment of  the  system  there  is  equal  probability  of  finding  either  molecule 
in  the  higher  energy  state. 

That  this  resonance  situation  applies  to  collisions  in  which  the 
collision  parameter  b is  so  large  that  there  is  no  appreciable  deflection 
of  the  trajectory  is  indicated  by  the  theoretical  interpretation  of  the 
observed  presstire  broadening  of  microwave  spectra.  Experimental 
values  of  the  collision  cross  sections  for  the  self  broadening  of 
spectral  lines  of  polar  molecules  give  collision  diameters  which  are 
three  to  four  times  larger  than  the  kinetic  theory  collision  diameters. 
That  is,  whereas  most  kinetic  theory  diameters  are  in  the  range  of 
three  to  five  .Angstroms,  the  microwave  diameters  usually  lie  in  the 
range  of  ten  to  twenty  .Angstroms. 

16 

The  theory  of  pressure  broadening  Indicates  that  the  largest 
cnniribution  to  the  cross  section  for  polar  molecuJes  comes  from 
resonant  collisions,  and  that  the  observed  ncucrowive  collision 
diameters  may  be  interpreted  approxdmately  as  the*,  critical  value 
of  the  collision  parameter  b for  which  the  probability  of  transition 


15.  D.  Bohjn,  loc.  cit.,p.  477. 

16.  For  a summary  of  the  microwave  theroy  and  experimental  data, 
and  a guide  to  the  literature,  see  W . Gordy,  W.  Smith  and 
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during  the  encounter  is  one  half 

Thus,  for  collisions  in  which  the  collision  parameter  is  less  than 
the  microwave  value  of  the  collision  diameter,  we  are  justified  in  con- 
.Sidering  both  final  states  of  the  system  as  equally  likely.  For  collision 
parameters  larger  than  this  value,  the  angle,  of  deflection  in  most 
encounters  is  negligible 
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APPENDIX  I 

ALMOST  ADIABATIC  TIME  DEPENDENT  PERTURBA TIONS* 


When  the  Hamiltonian  of  a system  changes  slowly  with  the  time, 
we  may  expect  to  be  able  to  approximate  solutions  of  Schrosdinger 's 
equation  by  means  of  stationary  state  energy  eigenfunctions  of  the 
instantaneous  Hamiltonian. 

Let  the  Hamiltonian  ot  a system  be  H (x,  t),  where  t is  the 
time,  and  x denotes  all  of  the  space  coordinates.  Let  v (x,  d) 
be  the  stationary  state  solutions  obtained  by  setting  t = 6 , a constant. 
The  ( X,  0 ) are  given  by 


H 


( 


/ 2- 


/4^/’ 


£^(6) 


(A-1) 


Ws  assume  that  the  (x,  6)  form  a complete  orthonormal  set  of 
functions  in  the  space  coordinates  at  every  time. 

If  the  total  wave  function  of  the  system  is  known  at  zero  time, 
wt  mav  at  later  times  write 


d& 


(A -2) 


The  expansion  coefficients  A^  (t)  depend  on  time,  and  the  form  of  the 
exponential  time  factor  has  been  chosen  for  convenience  in  subsequent 
differentiations . 


The  following  derivation  is  essentially  the  same  as  those 
found  inD.  Bohm,  loc.  cit.  , p,  49'»,  and  L.  I.  Schiff,  loc.  cit.  , 
p.  207. 
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We  may  obtain  the  expansion  coefficients  (t)  by  requiring 
that  satisfy  the  second  Schroedinger  equation 


. , ,-y  ^ 

I-! 


(A-3) 


This  gives 


<?■  .• 


i (e.  .., 


(A-4) 


Multiplying  by  t)  0 and  integrating  over  all  x, 

we  obtain. 


0 - (ic, 


(A -5) 


where  we  have  used  the  orthonormal  properties  of  the  v^  (x,  t). 

This  result  may  be  simplified  by  rewriting  the  matrix  elements 

(v  I V,).  From  Eq.  A-1 
m'  j 

H v-j  - Ej 

Differentiating  with  respect  to  the  time,  we  obtain 


M V,  ^ H - Ei  ^ tj 


IT: 


(A-6) 


♦ 


If  we  multiply  by  ^ j)  » integrate  over  all  x,  we  get 

immediately. 


( v^l  Hi  Vj)  i-  = 


^j) 


(A -7) 


Since  H is  Ilermetian, 

(ir^l  Hi  in)  - (■uv.l'vv) 
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Thus  Eq.  A -7  gives 


( 0/“^  i iTj) 


[z  j — 


(A  *8) 


and  Eq.  A-5  may  be  written 

Lw  - X .A-„ 

S #/V>7 

From  here  we  proceed  as  in  the  usual  method  of  the  variation 

of  constants.  If  we  assume  that  at  t = -CD  the  system  was  :'.n  the 

pure  state  v , we  may  approximate  the  A.  ft)  as 
o - - j ' ' 


(t)  - 


so  that 
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The  probability  that  at  time  t the  system  will  be  found  in  the 
m'th  state  is  now 
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